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and here are some of 
the many reasons why: 


OUTSTANDING 
ALL-ROUND PERFORMA''‘¢ 


DD MAZ NE. © Fast, reliable hypergolic ignition 
: with most storable and high- 


ROCKET FUEL FACTS 7 Smooth, safe starting and shut- 


down transients 
Excellent combustion stability 


High resistance to liquid-phase 
decomposition and vapor detona- 
tion in injectors, coolant tubes 


High specific impulse: 


855 80/20 Fluorine/LOX 
844 Fluorine 

810 LOX 

290 Perchloryl Fluoride 

286 Nitrogen Tetroxide 

282 98% Hydrogen Peroxide 
280 Chlorine Trifluoride 

276 IRFNA 

(Theoretical Isp, see—1,000/14.7 psia, 
optimum expansion, shifting equil.) 


Excellent efficiency, high actvall,, 


Low viscosity, good pumping and 
fluid-flow properties 


High heat capacity, good coolant 
characteristics 


g ty Density comparable to hydrocar- 


q bon fuels 


Outstanding storability 
Ease of handling 


Simplified equipment design and 
conversion from other fuels 


Extensive background of testing 
and operational experience 


e Rapidly expandable production 
to meet large-scale requirements 


DIMAZINE. 


unsym-Dimethylhydrazine-UDMH 
the storable high energy fuel 
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Revolutionary RCA Magnetic Video Tape Recorder to in 
Speed Navigation Training of Submariners 


Aboard the nuclear submarine Sea Dragon, the first 
undersea magnetic video tape recorder will record 
and store data on under-the-ice characteristics from 
externally installed TV cameras. Upon return to 
base the recorded information will be displayed for 
the benefit of undersea service trainees. 


The RCA undersea recorder is a marvel of compact 
design (dimensions 20” x 20” x 100”). It nestles in 
a torpedo rack, and represents a 60% space reduc- 


tion over existing video tape equipment. | 


Among the exclusive RCA developments are: the 

now famous “Tiros” satellite recorder; a radar sys- 
tem designed to take the first pictures of a nose cone 
re-entry vehicle; a unique tape cartridge completely 

adaptable to any size recorder. For literature de- 

scribing new RCA defense and commercial products 

developments, write Defense Electronic Products, 

Radio Corporation of America, Camden, N. J. 


Out of today’s defense needs...tomorrow's electronic advances 


The Most Trusted Name in Electronics | 
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Scope of ARS JOURNAL 

This Journal is devoted to the advance- 
ment of astronautics through the dissemina- 
tion of original papers disclosing new scientific 
knowledge and basic applications of such 
knowledge. The sciences of astronautics are 
understood here to embrace selected aspects 
of jet and rocket propulsion, spaceflight 
mechanics, high speed aerodynamics, flight 
guidance, space communications, atmospheric 
and outer space physics, materials and struc- 
tures, human engineering, overall system 
analysis, and possibly certain other scientific 
areas. The selection of papers to be printed 
will be governed by the pertinence of the topic 
to the field of astronautics, by the current or 
probable future significance of the research, 
and by the importance of distributing the in- 
formation to the members ot the Society and 
to the profession at large. 


Information for Authors 


Manuscripts must be as brief as the proper 
presentation of the ideas will allow. Ex- 
clusion of dispensable material and concise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references, 
1 page of abstract and 12 illustrations. 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acceptable 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes or Technical Comments. 
They may be devoted to new developments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are published within a few months of 
the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication, according to the above-sta 
requirements as to subject, scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as a 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be ahi spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ posi- 
tions and affiliations in a footnote on the 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the 
printer is essential. Greek letters and unusual 
symbols should be identified in the margin. If 
handwritten, distinguish between capital and 
lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to be given 
as follows. For journal articles: Authors 
first, then title, journal, volume, year, page 
numbers; for books: Authors first, then title, 
publisher, city, edition and page or chapter 
numbers. Line drawings must be clear and 
sharp to make clear engravings. Use black 
ink on white paper or tracing cloth. Lettering 
should be large enough to be legible after 
reduction. Photographs should be glossy 
prints, not matte or semi-matte. Each illus- 
tration must have a legend; legends should be 
listed in order on a separate sheet. 

Manuscripts must be accompanied by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. 

Preprints of papers presented at ARS 
meetings are automatically considered for 
publication. 
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issile Family...Scions of Space Tec 
re * 


Science and technology, especially as they relate to missile art, have advanced further in the last six years than in the pre- | 
ceding six centuries. Any review of the many milestones successfully attained since 1954 reveals an epic of hard work, 
inventiveness, accomplishment, and singleness of objective. This single objective—the achievement of operational weapon 
capability at the earliest possible date—is being realized. 

The Air Force missile family including Atlas, Thor, Titan, and Minuteman, has achieved progress beyond expectation in 
a program unmatched for magnitude and complexity. 

Space Technology Laboratories has had the responsibility since 1954 for the over-all systems engineering and technical 
direction of these programs. STL’s scientific and technical management capabilities have not only helped to hasten the day 
of operational capability for Air Force ballistic missiles, but have also been applied in carrying out related space probe and 
Satellite projects. 

Scientists and engineers with outstanding qualifications find unusual opportunities for their skills and disciplines at STL. 
Positions on STL's technical staff are now available for those who wish to add a new dimension to their careers. Resumes 
and inquiries are invited. 


SPACE TECHNOLOGY LABORATORIES, INC. 


P. 0. Box 95004, Los Angeles 45, California ¢ STL 
Florida offers immediate opportunities to outstanding 


Assistant Test Supervisors. Write: Mr. George S. Cherniak, 
STL, P. 0. Box 4277, Patrick Air Force Base, Florida 
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Aerospace Corporation has been 


brought into being to serve the United 
States government by concentrating 
the full resources of modern science 
and technology on rapidly achieving 
those advances in space systems indis- 
pensable to the national security. 

fhe corporation is non-profit, will 
share the findings of its research and 
laboratory experiments with all 
appropriate organizations involved in 
the government’s missile-space pro- 
gram, and is not organized for 
manufacturing purposes. 

he immediate responsibility of 
Acrospace Corporation is to aid the 
United States Air Force in bringing 
about the best possible ballistic missiles 
and military space systems on a con- 
tinving basis and within the shortest 
possible time. 

in addition, it may furnish the 
National Aeronautics and Space Ad- 
ministration and other governmental 
agencies appropriate services whenever 
its participation in space exploration 
and related activities is desired. 

'HE MISSION of Aerospace Cor- 
poration, according to the Secretary of 
the Air Force, encompasses “the field 
of ballistic missile and space programs. 
Within this complete area, it has the 
responsibility for advanced systems 
analysis, research and experimentation, 
and initial systems engineering. It will 
also exercise such general technical 
supervision of ballistic systems as is 
appropriate. In special cases, and with 
the consent of the Secretary of the Air 
Force, Aerospace Corporation may 
assume broader responsibility for an 
Air Force military system.” 

The new corporation also provides 
support to the Air Force in its effort to 
achieve maximum interchange of 
knowledge with other military services 
and among universities, research foun- 
dations, and the scientific community 
in general. 

THE FUNCTIONS of Aerospace 
Corporation in carrying out its 
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responsibilities for the Air Force’s 


missile-space programs include the 
performance of a wide range of 
scientific, technical, and administrative 
tasks. 

The corporation will conduct exten- 
sive laboratory and field activities 
aimed at advancing the state-of-the-art 
and will augment these research and 
development activities by coordination 
with industry, universities, labora- 
tories, and other agencies. It is in- 
tended that this combined effort will 
push forward the boundaries of tech- 
nology on a broad front to fulfill mili- 
tary and other national requirements. 

Aerospace Corporation will study 
the application of the advancing tech- 
nology to military weapons, support 
systems, and other systems serving the 
national need. These studies will 
culminate in preliminary design and in 
recommendations for development 
programs. 

Aerospace Corporation will then 
assist the Air Force or other appro- 
priate government agencies in es- 
tablishing space programs and in 
bringing the force of American indus- 
try to bear in carrying them out. Once 
development is initiated, Aero- 
space Corporation will assume 
responsibilities for the broad technical 
aspects of these new programs through 
their critical phases. 

THE FACILITIES of Aerospace 
Corporation include a research and de- 
velopment center located near the Los 
Angeles International Airport and 
within easy reach of several attractive 
residential communities. 

They constitute a modern adminis- 
trative, scientific, and engineering 
headquarters which house some of the 
world’s most advanced instrumentation 
and experimental apparatus. 

In addition, the operations of 
Aerospace Corporation will be directly 
supported around the globe by a vast 
array of resources created by the 
government over the past six years. 


2400 East El Segundo siiabininit El — California 


These will include: The Atlantic 
Missile Range in Florida; The Pacific 
Missile Range in California; The 
Rocket Engine Test Site at Edwards 
Air Force Base in California; and 
numerous other missile test facilities 
sponsored by the government in 
cooperation with private industry. 

THE PEOPLE who make up Aero- 
space Corporation have been selected 
from industry, universities, and 
government. They constitute a cross- 
section of highly-developed engineer- 
ing and scientific skills in the missile 
and space fields. 

Extensive recruiting will continue as 
the new corporation assumes more and 
more responsibility. From the outset, 
the corporation is built on a foun- 
dation of proven scientific competence, 
imagination, and objectivity. 

THE OPPORTUNITY awaiting 
those scientists and engineers who 
qualify to join Aerospace Corporation 
is equalled only by the magnitude of 
the corporation’s mission — magnitude 
mirrored by the highly advanced na- 
ture of the programs in which Aero- 
space Corporation is engaged. 

Typical systems projects include: 
advanced ballistic missiles; advanced 
military space boosters; recoverable 
boosters and satellites; space defense 
systems; early-warning satellites; 
reconnaissance satellites; communica- 
tions satellites; and manned satellite 
systems. 

Typical research programs concern: 
nuclear propulsion; astrodynamics; 
magnetohydrodynamics; inertial ele- 
ments; millimeter waves; hypersonics; 
combustion kinetics; and materials 
research. 

Those capable of contributing to 
state-of-the-art advances in these and 
related areas are invited to consider 
the advantages of becoming a part of 


the new Aerospace Corporation. Their 


resumes should be directed to: 
Mr. James M. Benning, P.O. Box 


95081-I, Los Angeles 45, California. 
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to meet the needs of three major projects 


MHD POWER GENERATOR—In cooperation with 10 leading 


joe sae utilities, AERL is working toward a prototype electrical power 


generating station using the revolutionary Magnetohydrodynamic 
technique. 


— 2, RECOVERABLE SATELLITES—Under contract with the U.S. Air 
ir Force, AERL is developing and flight-testing recoverable satellites 
based on the unique Avco Drag Brake principle. 


ICBM. DEFENSE—AERL's superior background in laboratory and 
Dateien: aT inflight studies of re-entry phenomena is being applied to the vital 
ae anise, 12 problem of the design of practical systems for defense against 

ICBM's, under several contracts with the Department of Defense. 


CONTINUING RESEARCH PROGRAM—These ticular, we are looking for: Aerodynamicists with 
expanding projects are results of AERL’s basic some background in high-temperature gas pheno- 
research effort. To provide general support for mena, Physicists and Physical Chemists interested in 
existing projects, and to supply new fundamental plasma and associated high-temperature gas 


knowledge for our further growth, we are seeking phenomena. ow d 
scientists and research-oriented engineers. In par- Mveo Sons a 
LABORATORY 


TO SPEED THESE PROJECTS, WE ARE SEEKING: 


Aerodynamicists— Communications & FOR FURTHER INFORMATION 
Theoretical & Experimental Telemetry Specialists ti 

Senior Rocket Combustion Engineers Physicists—Theoretica! & _ write to: 

Mechanical Engineers Experimental Mr. R. E. McDonald, 

Engineering Physicists Aeronautical Engineers 

Systems Engineers aay Missile and Satellite Designers Aveo-Everett 

Flight Test Specialists Physical Chemists “9 Research Laboratory, 

Electronic & Electrical Engineers Flight Control Specialists —_. Everett, Massachusetts. 
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Particle Impacts on the Melt Layer 
of an Ablating Body 


An analysis has been made of the effects of high velocity particle impacts on a melting body, such 
as might occur if atmospheric particles are encountered during atmospheric entry. The surface 
material considered in the analysis is a composite of glass and a plastic which evolves a large quan- 
tity of gas. An equation is developed relating the increase in erosion rate, as a result of particle 
impacts, to flight and material parameters. It was assumed in the derivation that the liquid melt 
layer protects the underlying solid body by decelerating the particles. The extent of the melt layer 
protection was found to be dependent on the particle size and melt layer thickness. One experi- 
mental check on the analysis was possible by inspection of the order of magnitude of the experi- 
mental liquid melt layer thickness on fired Micarta specimens. The computed melt layer thick- 
ness agreed with that observed on glass-reinforced Micarta 259-2 test specimens which were exposed 
in a rocket exhaust jet. Finally, it was assumed that the effect of the particle impacts was to trans- 
fer energy to the plastic, thus increasing the decomposition rate of the plastic filler in Micarta 259-2. 
It was possible to compute erosion rates, based on this analysis, which were in agreement with ex- 
perimental results. The analysis shows that the particle impacts influence the thermal energy 
tot body. 


VEHICLE entering a planetary atmosphere will be nute matter exists between 35,000 and 50,000 ft in the Martian 
subject to extreme heating. The extent of this heat- atmosphere; this matter is likely to be ice crystals. There 


E. W. UNGAR? 


ing, which results from dissipation of vehicle potential and is also the strong possibility of the presence of significant 
kinetic energy, is dependent on initial entry velocity, nature quantities of dust in the atmospheres of Venus and Mars. 
of the atmosphere, gravitational attraction of the planet and The Martian dust may be in the form of dust clouds. Thus, 
the weight to drag ratio of the vehicle. Entry into the at- in addition to the purely theemal consideration of atmospheric 
mospheres of Mercury and Mars will be less severe than entry entry, some thought must be given to the effects produced 
into Earth’s atmosphere. Entry into the atmosphere of by the impact of high velocity vehicles: upon atmospheric 
Venus is comparable with entry into Earth’s atmosphere; particles. These impacts are likely to occur during the 
whereas entry into the Jovian atmospheres will probably be aerodynamic heating phase of the entry mission. The 
much more severe. In order to relieve aerodynamic heat- effect of these impacts might be to increase the ablation rate 
ing, the leading edges of the vehicles used for missions re- to such an extent that the interior temperature of the vehicle 
quiring atmospheric entry will usually be blunted so that a would exceed the limit allowed for personnel or equipment. 
large amount of the vehicle energy will be dissipated by heat- Furthermore, the increased rate of ablation might occur at 
ing the gas behind a detached shock wave. However, there locations where only small amounts of material loss would 
still may be sufficient energy transferred to the body from otherwise be expected. In an extreme situation, the geometry 
the hot gas behind the detached shock wave to vaporize most of the vehicle might be altered sufficiently by particle im- 
materials. Thus, a thermal protection system is required pacts to affect the flight path. In the other extreme, if the 
so that the interior of the vehicle remains at a sufficiently low particles are composed of water, it is possible that the effect 
temperature to allow operation of equipment and the well- of the particles might be to cool the surface, resulting in a net 
being of passengers. Past experience with ballistic missiles decrease in the ablation rate. 
has shown that ablating materials are practical and reliable This paper presents the results of a study of the impact of 
for use in thermal protection systems. small particles on ablating bodies. The experimental ap- 
There are at present many uncertainties with regard to the proach and test results are briefly outlined. The impact 
exact nature of the various planetary atmospheres. Atmos- phenomenon is then analyzed, and the experimental data 
pherie moisture particles in the form of clouds and rain are correlated with the proposed theory. Finally, implications 
present in Earth’s atmosphere. There is evidence of moisture of the study with respect to atmospheric entry — are 
in the atmospheres of Venus and Mars (1). A layer of mi- briefly discussed. 
Presented at the ARS Structural Design of Space Vehicles General Research Procedure 7 
Conference, Santa Barbara, Calif., April 6-8, 1960. ; 

! The investigation reported herein was carried out for the The extent of the damage that might be caused by impact 
ad Ballistic Missile Agency under Contract DA-33-019-ORD- of a melting aerodynamic surface with suspended particles 
277 

? Principal Mechanical Engineer, Aeronautics and Thermo- during an atmospheric entry is influenced by - number of 
dynamics Division.:id¢Member ARS. variables. Among these variables are particle size, particle 
Numbers in parentheses indicate References at end of paper. density, impact velocity, surface material and surface condi- 
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tion. In addition to these variables, there are additional 
variables associated with the effects of the shock wave and 
shock layer gas on the particles, namely, particle deceleration 
and deflection, particle deformation and particle evaporation. 
An investigation of the effects of the shock layer gas on 
moisture particles has shown that, for a first-order approxi- 
mation, they may be neglected in the case of a vehicle entering 
Earth’s atmosphere (2). Therefore, the effects of the shock 
layer gas on the particles are not considered in this paper. 

Although perfect simulation of all pertinent variables is not 
feasible in a ground test facility, a reasonable approximation 
of the pertinent variables can be obtained in a rocket engine 
test facility. Rocket engines provide a stream of high 
velocity, high temperature gas into which particles may be 
introduced. High stagnation temperatures, of the order 
of 3000 to 6000 R, are easily obtained through the release 
of chemical energy within the combustion chamber. Al- 
though these temperatures may in some cases be lower than 
those at the atmospheric entry condition being simulated, 
they are high enough to melt most materials. Exhaust 
Mach numbers between 2 and 4 are obtainable at reasonable 
chamber pressures for sea-level operation. Although the 
Mach number in a rocket engine exhaust at sea level is gener- 
ally not as high as those experienced by a vehicle entering 
an atmosphere, it is sufficiently high to produce a flow, 
geometrically similar to that obtained in flight, about a 
properly located blunt body. The initial relative velocity 
between the particle and the shock layer gas will usually be 
less in the rocket engine facility than in flight, since rocket 
engine exhaust velocities are of the order of 7000 fps at sea 
level, resulting in particle velocities less than 7000 fps. En- 
ergy and momentum considerations, therefore, require par- 
ticles of greater mass than atmospheric particles. 

In summary, when the effect of the high velocity particles 
on a melting body is considered, a rocket engine test facility 
shows appreciable merit. Although impact velocities are 
not so great as those which might occur in flight, the particle 
energy or momentum can be simulated by adjusting the 
particle density. Furthermore, the stagnation temperatures 
and heat fluxes available in a rocket engine exhaust jet are 
high enough to maintain a reasonable test specimen thermal 
ablation rate. By controlling the size of the particles used 
in the experiment, it is possible to maintain the relationship 
between particle size and melt layer thickness required for 
simulation. Solid or liquid particle impacts might be ob- 
tained by varying combustion temperature or particle 
material. 

To gain some understanding of the physical phenomena 
associated with the impact of high velocity particles on an 
ablating body, experiments were conducted at the Army 
Ballistic Missile Agency. The facility used in these studies 
was the ABMA 3.5 AHT rocket engine. This engine burns 
a mixture of kerosene and hydrogen peroxide and is designed 
so that a slurry of additional kerosene and fine particle can 
be injected directly into the combustion chamber. The 
particles are accelerated by the hot exhaust gases and impact 
on test specimens. 

Fused aluminum oxide particles were used in a large 
number of the tests in the 3.5 AHT exhaust. Since these 
particles remain solid during passage through the engine, 
they were used to simulate the impacts of solid ice particles. 
Particles composed of a low-melting-point glass were used 
in the remainder of the tests. The glass particles melted in 
the exhaust jet and thus simulated the impact of liquid 
droplets. 

Micarta 259-2 is the material studied most extensively 
under this program, since it is considered to be typical of the 
thermal shield materials used on some atmospheric entry 
vehicles. Thermal ablation of Micarta in the 3.5 AHT ex- 
haust took place as a result of the melting of the glass lami- 
nates and the decomposing of the melamine resin filler. The 
melting glass resulted in the generation of a thin layer of 


flowing molten glass on the surface of the specimen, whereas 
the thermal decomposition of the plastic filler resulted in the 
blowing of gaseous decomposition products through this 
molten glass into the gaseous boundary layer. 


Experimental Results 


The exposed surfaces of the specimens were first studie:! 
under a binocular microscope using powers of magnification 
from 2 to 150X. Comparisons were made between speci- 
mens which were fired both with and without particles in the 
hot gas stream. Particles could not be identified by this 
type of examination. The solidified melt was observed to 
become more dense as the particle exposure became more 
severe. Some unmelted glass fibers could be seen in the 
melt layer on specimens exposed to particle impact; this wis 
not observed on control specimens fired without particles. 

Samples were scraped from the melted layers of repr: - 
sentative specimens and analyzed microscopically. Samplcs 
from the specimens impacted by aluminum oxide particl s 
showed aluminum oxide was present in the melt as discrete 
particles. The index of refraction of the melt was only 
slightly affected by the presence of the oxide. In addition, 
the size of the alumina particles found in the samples wei» 
measured, and it was found that the average size of these 
particles was less than the original size by 35 to 50 per cen‘. 
Cracks were also noticed in these alumina particles. The 
inferences to be drawn from these observations are that: 

1 The alumina did not melt in the rocket engine. 

2 Many of the particles fractured on impact. 

Thus, it is felt that solid particle impacts were adequately 
simulated by alumina in the 3.5 AHT facility. 

Samples taken from the melt layer of specimens whic! 
had been impacted by the low-melting-point glass par- 
ticles in the 3.5 AHT facility were examined. These ex- 
aminations showed that particle melting had occurred prior 
to impact on the specimens during the tests. It is to be 
recalled that liquid particle impacts were desired so as to 
provide simulation of liquid-water impact. The reason for 
believing that the glass particles impinged on the specimens 
while in the liquid state is that the microscopic examination 
indicated an intimate mixing of the glass from the Micarta 
and the glass particles. If the particles had struck the sur- 
face as solids, then the last particles to impact as the speci- 
men was removed from the test would be found as discrete 
particles. The melt layer on the specimens impacted by 
glass particles had indexes of refraction which were con- 
sistently less than either the index of refraction of the glass 
particles prior to the tests or the index of refraction of the 
melt on a specimen which was not exposed to particle impact. 
The decreased index of refraction is believed to be the result 
of some volatilization of the lead and boron constituents of 
the mixed glasses. 

In addition, observations of cross sections of the melt on 
several specimens were made. The cross sections were ob- 
tained by cutting the specimens with a diamond saw. Melt 
thicknesses were found to be generally of the order of 70 x. 
A charred layer was beneath the melt; its thickness was three 
to five times the melt thickness and generally thicker on 
specimens that were hit with particles. The particles 
might carry energy deeper into the body than would normal 
thermal conduction. Microscopic examination of samples 
scraped from the charred layer showed the presence of 
aluminum oxide particles in the charred layer of those speci- 
mens which had been impacted on by aluminum oxide par- 
ticles. 

Fig. 1 shows magnifications of the cross section through 
two typical test specimens. The melt layer on the specimen 
exposed to particle impact is considerably more nonuniform 
than the one on the specimen not exposed to particles. This 
is believed to occur during cooling. The melt is probabl; 
relatively uniform in thickness during the test; however 
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b 0.625 per cent concentration of 7.6-u glass particles 
Fig. 1 Cross sections through melt layer on specimens 


the melt layer on the specimen exposed to particles is prob- 
ably thicker than on the control specimen during exposure. 
Several unmelted glass fibers can be seen on the specimen 
which was exposed to the particle-laden gas stream. 

Fig. 2 shows the fractional increase in ablation rate due to 
particle impacts as a function of particle diameter for the 
Micarta test specimens hit by aluminum oxide particles. 
For each test the ablation rate obtained with particles m was 
normalized with respect to the ablation rate with no particles 
mo corresponding to the conditions of motor operation used 
for that test. Two different severities of test were utilized, 
yielding mp values differing by 50 per cent. The additional 
erosion resulting from particle impacts is seen to increase in 
a nonlinear manner with both increasing particle size and 
increasing particle concentration. Fig. 3 shows the particle 
impact velocity as a function of size calculated for conditions 
in the 3.5 AHT experiments. It shows that the absolute 
particle velocity decreases with increasing particle size. 
Therefore the total energy and momentum transport for a 
fixed weight concentration of particles decreases with in- 
creasing particle size. Thus, at first glance, an increase in 
the particle size with a fixed weight concentration of particles 
would not be expected to cause an increase in ablation rate. 
Fig. 2, however, shows that an increase in ablation rate was 
obtained with the larger particle sizes. This cannot be ex- 
plained solely on the basis of an increasing impaction eff- 
cieney with increasing particle size, because it has been shown 
that the impaction efficiency in these tests was close to 100 
per cent (2). Extremely large variations in impaction effi- 
ciency would be required to explain the results shown in 
Fig. 2. 

The cause of the nonlinear variation of the ablation rate 
with concentration was also not immediately apparent, since 
both particle energy and momentum transport vary in a 
linear manner with particle concentration. These experi- 
mental data are discussed further in the next section. 
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Fig. 2 Fractional increase in ablation rate for test specimens 
and aluminum oxide particles 


A melt protection mechanism is proposed to explain the 
results presented in the previous section. In some respects 
the proposed protection mechanism is similar to the pro- 
tection offered by thermal ablation. The mechanism is 
such that an increase in the ablation rate resulting from par- 
ticle impact results in a thicker melt layer. This thicker 
melt layer offers resistance to the particle and decelerates 
it, thus decreasing the impact energy and momentum of the 
particle as it strikes the solid underlying surface. A corre- 
lating equation based on this mechanism is developed. 

During the course of analysis based on the melt protection 
mechanism, it was found that the particle kinetic energy 
which is not dissipated in the liquid melt layer may be 
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accounted for if it is assumed that a fraction of this energy 
is used to decompose the melamine plastic filler in the Mi- 
carta. The decomposition of the plastic results in the 
evolution of a large quantity of gas. This gas has an effect 
on the thermal flux into the body. The energy required of 
the particles is only that necessary to start the melamine de- 
composition, because the decomposition itself is exothermic. 
The fraction of the impact energy which is transferred to 
the melamine plastic filler depends on the probability of a 
particle striking it rather than on the glass reinforcing 
fibers. This means of accounting for particle energy can be 
checked because the energy dissipated in the melt layer can 
be computed independently of any particular failure assump- 
tion. 

After passing through the melt layer, the energy of the 
particle is not sufficient to melt all of the glass fiber material 
eroded from the surface of the specimen. The microscopic 
observation disclosed the presence of unmelted glass fibers 
in the melt. The quantity of unmelted glass found in the 
melt region generally increased as the severity of particle 
erosion increased. This observation supports the theory 
that the particles break down or decompose the melamine. 
Additional support comes from the observation that the layer 
of charred plastic on fired test specimens was generally 
thicker when the specimen was hit by particles. The par- 
ticles might extend the decomposition of melamine to greater 
depths than would ordinarily be accounted for by heat 
transfer. 

The _ to be analyzed is characterized schematically 


in Fig. 
The ablation rate m is given by 


m = my {1] 


A quasi-static state is assumed; that is, the melt thickness 
is assumed not to be a function of time. The origin of the 


coordinate system is assumed to be on the solid surface and 
ablation. 


thus recedes with the surface durin 
quasi-static state the simplification 


For the 
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Fig. 4 Schematic model of the ablating system. Dimensions 
in the y direction are magnified 
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may be applied to the ablating system. Therefote, for the 
quasi-static state 


m = m1 
It is now assumed that m, may be divided into two com- 


ponents, a mass loss due to particle erosion and a mass loss 
due to thermal failure of material. Therefore 


m= My [2] 


On the basis of the assumed failure mechanism for the rein- 
forced plastic the following energy balance may be written 


Hfm, = F (KE) impact 


where F, the fraction of the particle impact kinetic energ 
transmitted to the plastic filler, may be taken as the volum 
fraction of plastic material in the composite reinforce 
plastic. Thus, F is essentially the probability of an impact- 
ing particle directly i impinging on the plastic as opposed to 
the glass laminates. 
The energy balance may be rewritten as ; 


1 


Rearrangement yields 
43 
m, = (F/2g.Hf)nw,Vs* 


The impact velocity at the solid surface must now be foun 
as a function of the melt layer thickness 6 and the initial par- 
ticle velocity V;. It is assumed that the particle passes 
through the melt layer before the cavity it forms as a wake 
has time to close. This last assumption is consistent with 
the order of magnitude of the melt thickness observed during 
the microscopic analysis of specimens. The particle velocity 


V is assumed to be normal to the surface of the body and not 
Newton’s 


to vary with distance from the stagnation point. 
law 


becomes for this case of cavity flow 


=» 


Algebraic manipulation yields 


dV _ V? 
dt 4 p,d, 
Equation [4] can now be written as 
2 
dV 3 pi V 
dy 4 p, ad, 
or as 
V 4 ppd, 
Integration of Equation [5] from V; to V; and from 4 = dto 
y = 0 yields 
= 


The value of V; will, in most cases, be taken as the free stream 
velocity of the particle relative to the body. The expression 
for the material loss due to particle impact may now be written 
as 


and the total ablation rate is rr 


bi 
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Expressions relating the melt thickness 6 and also ma to 
known parameters are still required. The thermal material 


loss my, will be considered first. Because there is an inter- 
action between the particle impacts and thermal effects, 
m,», is not equal to mp. It is now assumed that m,, is equal 
to the ablation rate mo, obtained when no particles are present 
in the free stream, minus the mass which would have been 
ablated by the energy used up in heating the gases, resulting 
from additional decomposition of the plastic, to the mean 
boundary layer temperature. Thus an interaction between 
particle and thermal effects is considered. The effect of these 
additional gases on the boundary layer thickness is not con- 
sidered. The interaction may be expressed approximately as 


— 0.6 (T. — Ty)e.fm,/He [9] 


where the 0.6 factor arises from Prandtl] number effects in 
th boundary layer. Substitution of Equation [7] into 
Equation [9], algebraic manipulation, and finally substitution 


int Equation [8] yield 


—3p,8/2p,d, 


mr = M 


[10] 


29.J Hf 
The only term now remaining to be evaluated in Equation 
[10| is the average liquid melt thickness 6. It will be assumed 
that all of the gas generated during the ablation process leaves 
the melt by flowing in the y direction, and all of the liquid 
and solid in the melt flows in the z direction. Following the 
example of Lees (3), the conservation of mass for the melt 

layer for axially symmetric flow is 
d 

2aro(1 —f)m = de ann dy 


where x is measured along the surface, and y is the coordinate 


normal to the surface. 
For an incompressible liquid melt + = = 


where A = 3 for a linear variation of velocity in the y direc- 
tion within the melt. 
Therefore continuity may be expressed as 
ro(1 — f)m = (dro/dx)Ap1uid [13] 
assuming that for small distances along the melt pz, u; and 
do are not functions of z. This assumption is not considered 
to be critical because large variations are not expected and 
because the effect of considering these quantities as vari- 
ables will tend to cancel when ratios are taken later in this 
discussion. 
For a spherical body 
= Rsin 0 


The function dro/dx may be expressed as 


dry _ dry do 
dx d0 dx 
but for a sphere 7 
= Rcosé 


which leads to 
dro/da = 


The equation of continuity may now be rewritten as _ 
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[14] 

The common assumption is now adopted that the shear 
stress is not a function of the melt thickness and would be the 


same if the melt were not present. This implies a linear 
velocity gradient in the y direction through the melt. Thus 


[15] 
[16] 


R tan (1 — f) m = Aprusd 


R tan 6 (1 — f) m = (1/2)prud = pitti 
105/50 
Substitution of Equation [16] into Equation [15] sila 
R tan 6m(1 — f) = 
of 


and solving for 6 yields 


UL 


~ 


[17 


It may be noted that for the condition of no particles in 
the free stream 


= V(1— f) 5omR tan 


tan — f)mo = pritirodo 


which may te solved for pxiizo to yield 
priito = m (1 — f) R tan 


which may now be substituted in Equation [17] to yield 


[19] 
An expression has now been obtained for the melt layer 
thickness which may be substituted into Equation [10] to 


yield the final result 


m mo E 0.6(7. | 


e 
1W,V 2p,d 


[18] 


5 = 


P 


PF 
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This equation may also be expressed in logarithmic form as 


"Vi — 0.6(7. — T;)fe,/H.] Fnw,V 2 
[21] 


Equation [21] implies that a semi-log plot of (m — mo)/ 
monwW,V 2 as a function of (m/mp)'/2/d, should correlate the 
data for specimens in the 3.5 AHT exhaust jet for one par- 
ticle density and for one thermal ablation rate. The negative 
of the slope of the resulting straight line should be the factor 
P,/P15o. Knowing this factor for one particle density and 
thermal] ablation rate should allow the determination of the 
factor for other conditions by adjustment of p, for other 
particle materials and 6o for other thermal ablation rates. 

Fig. 5 is the semi-log plot just described. All of the data 
plotted in Fig. 5 are for specimens which were impacted by 
alumina particles, p, 4 gm/ce, and the same severity of 
thermal exposure was used, i.e., the same mp. The impaction 
efficiency » is taken as unity. The correlation of the data 
shown in Fig. 5 is considered to be good. A slope of —29.1 
u~' is obtained, which yields a value of 77.7 u-gm/ce for the 
product pz5o. If pz is taken to be of the order of unity, 
then 49 becomes 77.7 u, which agrees with the order of magni- 
tude of the melt thickness observed on cross sections through 
specimens. As expected, it is less than the 400 » mentioned 
by Sutton (4) for Pyrex glass exposed in a rocket exhaust. 
The value of 59 for Micarta in the 3.5 AHT should be less 
than that for Pyrex in some other facility because the 3.5 
AHT gives relatively low thermal ablation rates and roughly 
half of the ablated Micarta enters the melt as a gas rather 
than as a liquid. 

The equation used for the correlation of all of the Micarta 
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data is obtained from rearrangement of Equation [21]. This 
equation is 


Hf(m — mo) 
{1 — 0.6 T ;)fce,/H.| Fyw,V;? 


Equation [22] may be checked for all of the Micarta speci- 
men tests conducted in the 3.5 AHT facility, including the 
tests wherein glass particles were used. The density of the 
glass particles is 3 gm/ce instead of the 4 gm/cc for the alu- 
mina particles, and the glass particles strike the specimen 
surface while in the liquid state. The value of pzé0 must be 
adjusted for severity of test utilized. 

Fig. 6 shows a plot of the right side of Equation [22] as a 
function of the inverse of the left side of the equation. The 
plot is presented on log-log coordinate paper because of the 
range of p,d,/pz50 which was tested. Had the left side not 
been inverted and had the plot been made on semi-log coordi- 
nates, a straight line would have been obtained which passes 
through unity when pz40/p,d, is zero. Instead of this hap- 
pening, however, the inverse plot of Ez shown in Fig. 6 ap- 
proaches unity asymptotically as p,d,/pzé0 gets very large. 
The plot implies that for a fixed impact velocity and flow 
rate of particles, the ablation rate approaches a constant 
value as the particle size is increased. This should be ex- 
pected because it indicates that the effect of the melt layer 
on the particle energy becomes very small when the particles 
are larger than some certain size. The variable Ez may be 
physically interpreted as an energy ratio, being roughly the 
ratio of particle energy utilized in bringing about failure of 
material to the particle energy immediately outside the melt 
layer. 

All of the data obtained for specimens in the 3.5 AHT fa- 
cility were plotted in Fig. 6. The correlation achieved with 
Equation [22] is considered to be good. The impaction 
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Fig. 5 Correlation of erosion data for specimens impacted by 
alumina particles at concentrations of 0.625 and 5 per cent and 
with average sizes of 16, 52 and 88 » 
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efficiency was taken as 100 per cent, and the value of V; wa: 
taken to be the particle velocity at the specimen locatior 
neglecting deceleration in the shock layer. 

The results presented in Fig. 6 show that the effects o: 
the liquid glass and the solid alumina particles on Micart: 
are of the same type. That is, it does not appear to mak: 
any difference whether the particles are liquid or solid whe: 
they impact on melting surfaces with high velocities. 


Implications of Results 


The analysis shows that vehicle velocity, melt thicknes:, 
ablation rate and particle size and density are all importani 
factors in the determination of impact damage. Therm:! 
ablation rates less than those encountered during entry int 
Earth’s atmosphere might be encountered during entry int 
the atmospheres of Venus and Mars. The lower therm: 
ablation rates, and possibly higher entry velocities, woul: 
tend to maximize, relative to an Earth entry, the damag: 
produced by particle impacts. 

As a part of the program described in (2), some calcula- 
tions of impact damage were carried out for ballistic missil:- 
trajectories. It is of value to consider the magnitude of the 
calculated damage, as a frame of reference when considering 
the damage which might be inflicted on a vehicle entering an 
atmosphere in which particles are present. 

Equation [22] was used to compute the magnitude of the 
increased erosion due to impacts with cloud elements of 
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Fig. 6 Ablation energy parameter as a function of the particle 
size ratio 
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body coated with Micarta 259-2 following a ballistic tra- 
jectory. The impaction efficiency was taken to be 100 per 
cent, and particle evaporation and acceleration during 
passage through the shock layer was neglected. The cooling 
effect due to the water particles acting as a heat sink was, in 
most cases, neglected. The calculations indicated an 8 per 
cent increase in the total mass ablated for the body following 
a severe re-entry trajectory as a result of flight through a 
nimbostratus cloud. The cloud altitude was assumed to 
extend from 2.3 to 10.9 km. The mean particle size was 
taken to be 13.5 uw and the concentration was taken to be 
5 10? particles per cc. At one point along a second 
trajectory which was studied, calculations indicated that 
the instantaneous ablation rate of the re-entry body may be 
increased by as much as a factor of 2 during flight through a 
cloud. 

in the case of vehicles entering Earth’s atmosphere, the 
particle impact problem will generally become more severe 
as the trajectories are changed to include longer periods of 
tine during which hypersonic velocities are maintained at 
altitudes less than 12 km. This would come about with 
ballistie missiles having higher weight to drag ratios and 
with hypersonic low altitude vehicles. Particle impact 
efiects also become more severe as smaller nose radii are 
utilized for re-entry vehicles. Some of the materials cur- 
rently being investigated for high heat flux applications 
do not form a melt layer and therefore might be more vul- 
nerable to particle impact damage than those which do 
melt during ablation. 

In the case of space vehicles entering other planetary 
atmospheres, the magnitude of the particle damage may be 
greater than for a ballistic missile entering Earth’s atmos- 
phere, in a cloudy location, if the particle concentrations en- 
countered are of the same order of magnitude. A large 
amount of damage can result from high initial entry velocities 
and an atmosphere of lower density than Earth’s. Because 
of these two factors, the vehicle entering from space might 
be traveling at higher velocities at the altitudes where par- 
ticles might be encountered than velocities realized with 
ballistic missiles in Earth’s atmosphere. The lower atmos- 
pheric density would reduce aerodynamic heating, and, as a 
result, the protective melt layer might not be as deep as dur- 
ing a comparable Earth entry. 
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Nomenclature 


specific heat of the gases evolved from decomposition 


Cg = 
of the plastic 
dp = particle diameter 
Ep = ablation energy parameter defined by Equation [22] 
F = fraction of particle impact kinetic energy transmitted 
to the plastic filler 
Fy = force on particle 
f = mass fraction of plastic filler in the undamaged 
material 
Je = consistent conversion factor 
H = energy necessary to start the decomposition of the 
plastic filler 
i, = effective heat of ablation 
J = conversion factor to energy units 
Mp = particle mass 
m = total ablation rate 
mi = rate at which material crosses the liquid-solid interface 
Mp = portion of ablation rate due to parti: le impact 
Mstored = Yate at which the mass contained in the melt layer 
changes 
Mer = portion of ablation rate due to heat transfer 
ro = distance measured radially from axis of body 
l = time 
= free stream velocity 
u = velocity in z direction 
V = particle velocity 
Wp = mass flow rate of particles q 
= impaction efficiency 
0 = angle measured from stagnation point 
A = constant sera by Equation [12] 
Subscripts 
e = edge of gas boundary layer 7 
f = solid-liquid interface 
g = gas 
2 = gas-liquid interface 
L = liquid 
0 = no particle impact condition 
p = particle 
6 
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HE KNOWLEDGE of the ablation process around 

the surface of a blunt body is important for the optimum 
design of satellites, ICBM’s and any other bodies that must 
re-enter the atmosphere at hypersonic velocities. In this 
ablation process, the aerodynamic heat energy produced by 
friction is absorbed by the mechanisms of melting and 
vaporization of the surface material, mass injection into the 
air boundary layer from the vaporization of the material, 
and the radiation from the high temperature of ablation ma- 
terials that have good emissivities. An additional mecha- 
nism of heat absorption is present when the ablating materials 
react chemically with the air boundary layer. 

The aerodynamic heat transfer to the surface of blunt bodies 
that fly at hypersonic velocities inside the atmosphere has 
been solved in (2 to 5),3 which consider the case of laminar 
and turbulent heating without mass injection into the air 
boundary layer. The effect of mass injection on the laminar 
heating rates can be obtained from the solution of the corre- 
sponding differential equations (6,7). For the turbulent 
case, the theory is not so well developed as for the laminar 
case. However, the available theoretical and experimental 
results (8,9) indicate that the mass injection effect on turbu- 
lent heat transfer is less important than the corresponding 
laminar case for injection materials that have a molecular 
weight equal to or heavier than air. For materials that are 
lighter than air, the mass injection effect seems to be about 
the same for both laminar and turbulent boundary layers. 

The melting of the body surface creates a liquid layer under- 
neath the air boundary layer. The high aerodynamic fric- 
tion forces on the surface of such bodies—combined with 
minimum weight and high nose cone performance require- 
ments—call for an ablating material that has a very high 
viscosity and low thermal conductivity. High viscous 
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* The energy equation (see Eq. [44]) indicates that for a given 
surface temperature, high values of hers require high values of the 
fraction of vaporized material f. This fraction is a strong func- 
tion of viscosity. A material like Pyrex—which has material 
properties like quartz except for a lower viscosity—yields 
negligible values for f as compared with those for quartz even 
for the case of turbulent flow conditions. 


Ablation of Glassy Material 


The steady-state equations of motion for a thin layer of an incompressible glassy material on os 
surface of an ablating and radiating blunt body are reduced to a first-order ordinary differentia! 
equation which is integrated numerically. This solution is coupled with the solution of the ai 
boundary layer for both laminar and turbulent heat transfer with or without mass vaporization 0! 
the ablating material. The distribution of the effective energy of ablation around the body is thu 
obtained for a cone cylinder with a hemispherical] cap that re-enters the atmosphere at hypersoni: 
flight speeds, and has quartz as the ablating material. It is found that the ablation process fron: 
turbulent heating is more efficient than from the laminar case because of increased vaporization 
This solution of the equations of motion at the stagnation point has been verified by arc wind tunne! 
experiments. The present state of development of the are wind tunnel does not permit its use fo: 
experimental investigations of ablation around blunt bodies under turbulent heating. one 
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forces will prevent the molten layer from leaving the surfac« 
before its energy of vaporization has been utilized for the re- 
duction of heat energy into the body. The low conductivit; 
will of course retard the heat flow into the body, and thus 
yield a more efficient ablation process. The equations of 
motion for a thin liquid layer on the surface of a body are 
in general those of the incompressible laminar boundary 
layer, with the appropriate boundary conditions. The 
requirement of a high viscosity leads to further important 
mathematical simplifications of the Navier-Stokes equation 
for a thin layer (10,11) which then yield approximate solu- 
tions that are extremely important for practical applications, 


Equations for a Highly Viscous Liquid Layer 


The equations of motion for a thin nonradiating liquid 


layer on the surface of the body are (12): ote dl 


Continuity 
ra) 
(ur) + — (vr) = 0 (1 
oy 
Momentum 
Du dp Pa) ( 2) 


Energy 


oT. a/,ar dp 
tee 


Viscosity law 
= [4] 


In these equations, the coordinate system is the same as in 
the gas boundary layer. Because y = 0 at the wall or gas- 
liquid interface and is positive in the gas layer, y is negative 
in the liquid layer. These equations can be solved by the 
usual numerical procedure of laminar boundary layer theory 
{as was done by Sutton who considered the liquid layer of 
Pyrex (13)]. However, this numerical procedure is too 
cumbersome and time-consuming for practical applications, 
and any approximate method of integration and appropriate 
simplification of these equations have great practical im- 
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portance. Von Kérman devised such an approximate method 
of solution by integrating the equations of motion in the y 
direction (12). Also, for a fluid with a very large viscosity, 
the inertia forees become considerably smaller as compared 
with the viscous forces, because the former vary with the 
velocity squared while the latter are linear with the flow 
velocity. The magnitude of the Reynolds number becomes 
of order unity, and the particle inertia forces may be dropped 
from the momentum equation (12,13). This condition is 
approximated by liquid glassy materials that do not have a 
well-defined melting temperature (10). An experimental 
verification of this approximate method of solution has been 
obixined at the stagnation point of ablating models in the 
arc wind tunnel (1), where it was found that the maximum 
dis-repancy between theory and experiment is within 8 
per cent. 

‘hus, the equations of motion for a liquid, highly viscous 
ra(ating material that are to be integrated in the y direction 
are: 


Continuity 


re) 
az iy (vr) = 0 
Momentum 


Energy 


oT aT d 
T 
(ud? +097) zt oy tr) tu dy 


Viscosity Law 
u exp (AT-) [8] 


The term pG@ sin @ in Equation [6] represents the component 
of the missile inertia force on the liquid particle, and G is the 
number of g’s acting along the missile longitudinal axis. It 
should be noted also that the energy Equation [7] includes 
the radiation heat flux 7. The viscosity for a glassy material 
(Eq. [8]) can be expressed as an exponential in AT! (A = 
constant), which indicates that 4 —~ © as the temperature 
drops in a direction from the gas-liquid interface toward the 
‘interior of the body. 
The integration of the continuity equation [5] yields 


10 y 

or, using the conditions at the interface where ¥ = v,» and 
y= 0 
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Equation [10] is the expression for the fraction of ablated 
material v,/v_.. that is vaporized into the gas boundary 
layer. It will be seen that this parameter together with the 
vaporization velocity v, and the wall temperature 7’, fix 
the effective enthalpy of ablation, which is defined as the 
energy that the ablated mass rate must absorb in order to 
balance the aerodynamic heat rate without mass injection 
effect go, i.e., pv-ohers = go. The lower limit of integration 
is taken as — in the foregoing equations because the ve- 
locity u approaches 0 at the solid surface of the ablating body. 
This velocity u is obtained from the momentum equation 
[6], which yields 


vy a 
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where Tv.(z) is the gas shear at the gas-liquid interface, 
and can be evaluated approximately from the corresponding 
heat rate through the Reynolds analogy. The lower limit 
of integration is taken as —@ because the viscosity » ap- 
proaches in the interior of the body. Ina limited tempera- 
ture range, the exponential form of the viscosity law can also 
be represented empirically as 4«7'™, where n is a con- 
stant that depends on the temperature range, i.e. ,n «7, 
(10,14). Therefore, the following expression may be obtained 


B= pA T/T.) [12] 


where n is large compared with unity. For example, for 
Pyrex, nis approximately 6 at 7 = 4000 F, and for quartz 
nis approximately 20. Since in the energy equation [7] the 
predominant terms involve the temperature gradient in the 
y direction, it is reasonable to assume that at any point around 
the surface of the ablating body the viscosity variation through 
the liquid layer will be as at the stagnation point, where 
u«e¥’® (10). However, the full energy equation [7] is now 
used to relate the ablation velocity v_.. to the liquid layer 
thickness 6. Thus 


[13] 


where 6 = 6(z). This expression indicates also an increase in 
viscosity in the liquid layer as y ~ —@. Combining Equa- 
tions [12 and 13] 


= T, [14] 
From which 
noT/dOy n(OT/OY)o Que 


Thus, the liquid layer thickness 6 becomes a boundary condi- 
tion. It is important to note that the heat rate qu is evalu- 
ated in the liquid side of the gas-liquid interface. From 
Equations [11 and 13] 


P,6? 
uly) = ult ¢ 1) [16] 


Mw 


which when introduced in Equation [10] yields 


Vv Tw? Ps 
=1- — — — 2— l 


The wall temperature is related to the ablation velocity 
through th uation [7 
gh the energy equa = 
0 oT oT ~ 
ut 


1 1 0 
— — +— (up: + dy [7a] 


ply 
or 
where 
gi = — 


of 1 fo ar/a 


0 
= + wey”) dy [22] 


The integration of Equation [7a] includes the heat conduction 
process in the interior of the body because there u approaches 
0, v = v_.. and all the terms except the second term in the left- 


: 
3 
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hand side and the first term in the right-hand side of the 
equation become, then, negligibly small. The conditions 
of negligible temperature and thermal heat flux have been 
used as boundary conditions at the interior of the body. 
The heat rate q; as defined by Equation [19] is the aero- 
dynamic heat rate with mass injection effect (¥qo) reduced 
by the energy required to vaporize the ablating material 
that is injected into the gas boundary layer (pzv.h.). The 
function y is unity for zero mass injection effect; otherwise 
it is reduced from unity by a term that depends on the nature 
of the heating rates (laminar or turbulent). Equation [18) 
is the steady-state heat balance equation. Since the vapori- 
zation velocity v» in Equation [17] is determined from the 
diffusion process in the gas boundary layer (Eqs. [31,32]), 
Equations [17 and 18] are the solution to the steady-state 
ablation problem, for they yield the ablation velocity v_. 
and wall temperature 7... Thus from Equation [18] 


Vi» = CE — — |/pepT w [18a] 


where ¢v_,, is obtained from Equations [14, 16 and 20] 


where conditions n=! << 1 have been used. The expression 
for Eo is 


—1 
—@ Tw 


which, with the aid of Reynolds analogy Shp: 7 
Ty = /*/(H, = H.) [23] 

and Equation [18] becomes : 


The condition of a positive tangential velocity at the gas 
liquid interface requires P.6/r. <1. It is thus seen that & 
has the order of magnitude of uwz/u. which is negligible com- 
pared with unity and may be neglected in Equation [18]. 
Except at sharp corners in the body contour where there are 
large Prandtl-Meyer expansions, £: is much less than 1. 
This term is usually small because it is produced mainly by 
the energy transfer in the liquid layer in a direction parallel 
to the surface of the body, and this is small compared with 
the energy transfer normal to the surface of the body. 

Equation [17] is a first-order differential equation that 
requires one set of boundary conditions, and can be solved 
more conveniently by using a stepwise numerical procedure 
that starts at the stagnation point of the body, where by 
symmetry 6. = Tv: = 0. This general equation [17] can 
be specialized to the stagnation point of a blunt body and 
to a semi-infinite cylinder as follows. 


Stagnation Point 


At the stagnation point region the conditions r = 2, 
Tw = Twrkt, Pr = bz » = 0 can be used in 
Equations [17 and 21]. The result is 


2 2 
8 


2 wad” SP... 


oN Twr et 


Cylindrical Body 


At the cylindrical section of an ablating body that is long 
enough to be unaffected by the ablation history ahead of it, 
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Equation [17] yields 


2 3 
CYL 


Vi« dz My 
( Pb Tw2d ) 


Effect of Radiation Mean Free Path on Ablation 


The general treatment of radiation effects on ablation 
is given in (15 and 16). It is shown there that the radiation 
mean free path has a strong effect on ablation, as can be 
seen from Equations [15 and 18], which yield 

pcpv_.n(1 é) qw 


When the radiation mean free path \ is small compared wit! 


T. Tw ns \2 tw n 


| 


— PA/re 1 — 1 — [ods — 


the liquid layer thickness 6, i.e., \/6 < 1, the liquid layer 
absorbs the radiation flux. In this case, a heat balance 
across a liquid sublayer having a thickness A at the gas-liquid 
interface yields gv — 0 = qi — eoT~‘. This is so because in 
the liquid side of this sublayer there are two equal radiation 
heat fluxes but in opposite directions. In the gas side of the 
sublayer, because the gas is transparent to radiation, there 
is only one radiation flux towards the gas boundary layer. 
Thus, Equation [28] becomes 


= k/pepnv-o(1 + &) (29 | 


It must be pointed out that in this case Equation [14] does 
not hold in this sublayer at the gas-liquid interface, and must 
be modified as indicated in (16). This is a local effect, how- 
ever, for Equation [14] is a good approximation for the 
temperature distribution in the liquid layer and interior of the 
body even for this case. On the other hand, when \/é6 > 1, 
the liquid layer is transparent to radiation. The conditions 
at the gas and liquid side of this sublayer are the same. 
Hence, the heat balance yields gq. — eoTy4 = qi — eoT.' 
Or Go = qi. The liquid layer thickness becomes then from 
Equation [28] 


(24) k ( 


Thus, the radiation mean free path affects the liquid layer 
thickness, and therefore, the fraction of vaporized material. 
Diffusion Equation and Vapor Pressure Law 


From the diffusion equation in the gas boundary layer in 
the absence of chemical reactions and neglecting pressure 
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gradient effects, is obtained (7) 


_ ¥qo (2 


Maie 


where the vapor pressure p, of the vaporized material must, 
of course, be known as a function of the wall temperature 
and the aerodynamic pressure. (For glassy materials where 
silica sand SiO, is the principal ingredient, see Appendix.) 
The parameter, Dz,/k in Equation [31] is obtained from 
the Lewis number, which for glassy materials can be as- 
sumed to be unity. 


Mass Injection Effect on Heat 


The injection of the vaporized ablated material into the air 
boundary layer reduces the heat transfer rate and the gas 
sheir at the gas-liquid interface. These two parameters 
can be related approximately by using the Reynolds analogy. 
For a given injection material and injection rate, the reduc- 
tion in the heat transfer rate depends in general on the type 
of flow in the boundary layer. 


ve 
Laminar Flow -s 


An investigation was made of the effect of air mass in- 
jection on the heat transfer from a partially dissociated air 
laminar boundary layer in thermodynamic equilibrium. As 
discussed in (3), the respective momentum and energy 
equations are 


(Ufnn)n + ffnn + B(p./p — fr?) = [33] 


+ fon + — fonln = 0 [84] 
The local pressure gradient and Mach number effects on 
heat transfer are given by 6 and u,?/H,, respectively. 
The mass injection effect is introduced as a boundary 
condition at the wall to the foregoing system of differential 
equations. From the results of (3) 


where the conditions of local similarity fe = 0 and negligible 

slip (fnw ~ 0) have been used. This latter condition is possible 
because the velocity u at the gas-liquid interface is very small 

compared with the velocity at the edge of the boundary 

are Introducing the blowing parameter B which is de- 
ned as 


2F 


vi 


= 
qo, V2 | Jnw s = [37] 


where (3) 


= 0.76 (Prptw) H, (1 GJws) Xx 

due 1/2 hos 


The combination of Equations [35 to 38] for a Lewis number 
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Effect of air injection on heating rates from a partially 
dissociated laminar boundary layer 


Fig. 1 
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of unity yield = & 4 


= [1.068/(1 + 0.096)VB)} [39] 


where 


= 


The ratio of heat transfer rate with mass injection to the 


heat transfer rate without injection ¥ as obtained from these 
numerical solutions can be correlated by the parameter shown 
in Fig.1. Itisseen that 


= 1 — 0.580"-4 + 0.096V8)/1.068] [40] 


This correlation holds for the injection of air into the disso- 
ciated air boundary layer. It shows a very weak dependence 
onl. The power of the Prandtl number is the result of two 


effects: 
1 The quantity 1 — y varies as o°-*. — 
2 The heat rate go varies as o~-°, a. 


The effect of pressure gradient and dissipation terms on this 
correlation was found to be small. There is also a somewhat 
weak effect of the molecular weight as is shown by the results 
of (7). A correlation of the results of (7) indicates that the 
blowing parameter B should be multiplied by A/a; 1/4, where 
Ma; is the ratio of the molecular weight of air to that of 
the injected material. Thus for laminar flow the mass in- 
jection correlation becomes 


= 1 — 0.5824 8M + 0.096V 8)/1.068 [40a 


Turbulent Flow 


The effect of mass injection on the turbulent heat transfer 
rates can be obtained from the experimental results of (8 
and 9). Reference (8) shows data for local heat transfer 
rates for the case of air injection into a turbulent boundary 
layer on a flat plate with an external Mach number of 2.7. 
Reference (9) presents data for the total skin friction coeffi- 
cient for the case of injection of helium, air and freon-12 
(CCI.F,) into a turbulent boundary layer on a cone with sub- 
sonic and supersonic flow at the edge of the boundary layer. 
This data indicates a weak Mach number effect on skin fric- 
tion for small injection rates. As the mass injection rates 
increase, the reduction in skin friction is larger at subsonic 
speeds. Using this latter data for the case of supersonic flow, 
the following approximate correlation is obtained for the total 
skin friction coefficient Cr Pia 


Cr ) 
fad’ Cr, 


[41] 


or 
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Fig. 2 Effect of mass injection on the turbulent skin friction 
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coefficient for a cone, M,. = 3.21 


a=4 when Ma <1 
a = when Ma; > 1 


This correlation is shown in Fig. 2. The relation between 
the local (C;) and total (Cr) skin friction coefficients is given 


where A is the wetted area of the body. Since at a given 
body location, the ratio C;/Cyo is independent of body posi- 
tion, Cr/Cr, is equal to C;/C;,. From the results of (5), 
Cr Hence, Cr, = C;.A, where 


Theoretical and experimental data indicate that both the 
skin friction and heating rates for turbulent boundary layers 
decrease with mass injection. These two parameters can 
be related approximately by extending the Reynolds analogy 
to the case of mass injection, in particular for the case Ma; < 


1 where the magnitude of this effect is smaller (9). Hence, 
Equation [41] can be written as 
Cy 2F 
= = Mat i 41 
OTPy)wo ~ Cr 
_or with Equation [36] 
5 


For « = 0.72, the parameter Ao’/* equals § for a cone and 
unity for a cylinder. The turbulent data of (9) show a 
large reduction in heat transfer rates for light injection ma- 
terials, i.e., materials that have a molecular weight lower than 
that of air. Equation [41] indicates this fact through the 
value of a = 4. For materials that are heavier than air, 
the reduction in heating rates is less important, as is indi- 
cated by a value of a = 4. 

It is also noted that for air injection (o = 0.72), the nu- 
merical coefficient in the laminar equation [40] becomes three 
times larger than the corresponding numerical coefficient in 
the turbulent equation [41b]. This indicates that for in- 
jection materials like air, the effect of mass injection on heat 
rates is less significant for turbulent than for laminar flow. 
Equations [40a and 41b] indicate, respectively, the effect of 
mass injection on the gas shear and heat rates for the condi- 
tions of laminar and turbulent boundary layers. 
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Determination of Ablation Parameters 


For a given blunt body geometry and flight conditions, the 
pressure distribution around the body can be determined 
independently of the ablation equations. The aerodynamic 
heating rates and local shear forces without mass injection 
effects for both laminar and turbulent boundary layers then 
become proportional only to the wall temperature of the 
ablating body. The ablation calculation is started at tie 
stagnation point where Equations [24 and 25] are uscd 
instead of Equations [17 and 20]. The use of Equation 
[17] around the surface of the ablating body requires a ste»- 
wise numerical procedure. By assuming a wall temperatui ., 
it is possible to calculate all the terms in Equation [24] «t 
the stagnation point, and in Equation [17] around the bo: y 
through the successive use of Equation [31] for ve, Equaticn 
[18] for v_., and [29] or [30] for 6. The viscosity is obtain: d 
from the corresponding equation [8], and the gas shear «t 
the interface from [23]. Note that at the stagnation poiit 
where the shear and velocity u, are proportional to x, Tw: is 
then proportional to the velocity gradient u... For a New- 


tonian pressure distribution 


An iteration procedure on wall temperature is then utiliz d 


— until either Equation [24] or [17] is satisfied. 


Effective Energy of Ablation 
The parameter that indicates the aerodynamic het 


energy that is absorbed per unit mass of ablating material 
is defined by 


hett = [42] 


where go is the aerodynamic heat rate per unit area without 
mass injection effect to the ablating body. This parameter 
heer is a function of both location on the body contour and time 
along a flight trajectory. An integrated energy of ablation 
for the ablation process during re-entry can be defined as 


ty ty 
hett = wat/ pv_ dt [43] 


The numerator is the total local heat input during the time 
of ablation, and the denominator is the total local ablated 
mass. This parameter herr is only a function of position on 
the surface of the body. From the energy equation [18], 
her: can be written as 


hes = {(1 + + + K(H, — X 
(1 — [44] 


where K is the coefficient of the blowing parameter in [40a] 
or [41b]. Equation [44] can also be expressed as 

heat = (1 + Te hy + (1 — het haw [45] 
Equation [44] shows the rather powerful effect of radiation 
on hers as q,/qo ~ 1, which is the limiting condition when 
all the aerodynamic heat is radiated away from the body. It 
should be pointed out, however, that this condition also 
gives a somewhat higher surface temperature, which yields 
higher g, and lower go rates. A higher surface temperature 
means a thicker insulating material for a given value of the 
temperature behind the heat shield material (17). Thus the 
effective energy of ablation as well as the surface temperature 
are the important parameters for the determination of the 
heat shield material along a given flight trajectory. 

Equation [45] indicates the breakdown or contributions 


to the effective energy of ablation by the mechanisms of hest« 


conduction, vaporization, mass injection and radiation, respec- 
tively. From this equation (or Eq. [18]) it is also seea 
that for steady-state ablation h.1: is always positive. 
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Applications 

The foregoing analysis is now illustrated by using a body 
shape of the Thor-Able type and typical trajectories for hyper- 
sonic flight (18). The ablating material is chosen to be quartz. 
The properties of this material are as given in the Appendix. 
A Newtonian pressure distribution is assumed for the hemi- 
spherical cap and cone, and a Prandtl-Meyer expansion is 
assumed to take place at the junction between the cone and 
cylindrical section of the body. For the sake of simplicity, 
the pressure on the cylinder is assumed to be constant. The 
heat rates are calculated by assuming that the air is in ther- 
modynamic equilibrium. It is, of course, necessary to know 
the conditions for transition from laminar to turbulent 


2.0 
STAGNATION PT 
----—SONIC PT 
| 
u | 
| | 
| 
12 > 
@ | | 
2 
x 
| \ 
o / \ 
\ 
7 \ 
7 
\ 
\ 
Oo 8 16 24 32 40 
TIME, SECONDS 
Fig. 3a Heating rates at stagnation and sonic pcints for typical 
5000-nautical mile trajectory. Ry = Roy./Ry = 1.5, ac = 
30 deg 
20 
15 
10 
= 
4 \ 
\ 
x = 
ONE——>| 
| 
fe) 0.4 0.8 1.2 1.6 2.0 
S/ Royt 
Fig. 3b Total aerodynamic heat input per unit area around body. 
Typical 5000-nautical mile trajectcry. Ry = 1% ft, Roy, =1.5, 
ac = 30 deg 
SEPTEMBER 1960 


= 


boundary layers. Based on the shock tube data of (19), transi- 
tion is assumed to take place at Re = 250. 

Fig. 3a indicates the aerodynamic heating to the ablating 
surface, without mass injection effect, during re-entry at the 
stagnation point and sonic point on the spherical cap of a 
typical nose cone that has a range of 5000 nautical miles. 
The sonic point shows the effect of transition from laminar 
to turbulent flow after about 14 sec of flight time in the 
atmosphere. It is also seen that the peak turbulent heating 
rate at the sonic point is about two times larger than the 
corresponding laminar rate at the stagnation point. The 
total heat input per unit area Jf qodt distribution around the 
body during the ablation process is as shown in Fig. 3b. 

The ablation velocity time history at the stagnation and 
sonic points of the above nose cone is indicated in Fig. 4a. 
This is the velocity with which the surface of the body moves 
toward the interior of the body. It is noted that the peak 
values of the ablation velocity at these two points differ by 
about only 15 per cent, even though the corresponding laminar 
and turbulent heating rates differ by a factor of 2.0. The 
reason for this is that there is more vaporization in the abla- 
tion process with turbulent heating than in the laminar case. 
The ablated thickness J y_.,dt along the body contour for the 
typical 5000-nautical miles nose cone is shown in Fig. 4b. 
The effect of emissivity on ablated thickness is also illustrated 
for the transparent case, i.e., \/6 > 1. It is seen that the 
maximum ablated thickness takes place near the sonic point 
on the spherical nose, where the turbulent heating rates 
reach a maximum along the body contour. The maximum 
change in body radius takes place when there is no radiation, 
and is about 7 per cent of the original radius. The loss of 
ablated material from the surface of the cone is less than at 
the stagnation point because of the increased vaporization, 
since heating rates are comparable (Fig. 3b). The ablation 
thickness on the cylinder is much smaller than at the cone, 
and is caused by the sharp drop in the heating rates at the 
cone-cylinder junction which is produced by the Prandtl- 
Meyer expansion in that region. In actual flight, as the pres- 
sure and, consequently, the heating rates, decay along the sur- 
face of the cylinder, the ablation process will eventually stop. 
Under these conditions, the assumption of steady-state abla- 
tion breaks down. The unsteady effects in the continuity and 
momentum equations for the liquid layer are negligible be- 
cause the density is constant and the inertia forces remain 
small compared with the viscous forces. However, the energy 
equation must be modified by retaining the term 07’/Ot in the 
left-hand side of Equation [7]. Also, these equations of mo- 
tion will not yield solutions for the combination of large 
inertia (G’s) and low aerodynamic shear forces. 


The integrated effective energy of ablation (her) distribu- 
tion around the body is as shown in Fig. 5. It indicates that 
for a given emissivity, the value of het: is considerably higher 
for turbulent flow conditions on the cone than for the laminar 
values at and near the stagnation point. As the emissivity 
increases, herr also increases because of the increase of energy 
that is radiated away from the surface, which produces mainly 
a decrease in the ablation velocity. The effect of both 
emissivity and radiation mean free path on the energy of 
ablation is as illustrated in Fig. 6. It shows that for a given 
emissivity, herr is higher for the case when the liquid layer 
is transparent to radiation. Fig. 7 shows that increased 
radiation tends to cool the ablating surface and that the radi- 
ation mean free path has a negligible effect on surface tem- 
perature. It should be mentioned that the viscosity data 
for quartz at high temperatures are not well known. How- 
ever, since the viscosity is an exponential function in the 
inverse value of the temperature, a plot of In uw vs. T-! is a 
straight line. This facilitates the extrapolation of viscosity 
data to higher temperatures. It was also found that if the 
viscosity at the highest temperature of interest (i.e., 6500 R) 
is varied arbitrarily by a factor of +3.0 from the extrapo- 
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Fig. 4b Ablated thickness around body for quartz. Typical 
5000-nautical mile trajectory. Ry = 1% ft, Reyt/Ry = 1.5, 
a. = 30 deg 


lated value (see Appendix), the total ablated thickness 
varies then by a corresponding +20 per cent. Fig. 8 indi- 
cates the contribution of the several mechanisms of heat 
absorption to the effective energy of ablation for the case 
where the emissivity is 0.5. It is noted that as the gas bound- 
ary layer becomes turbulent near the stagnation point, the 
fraction of vaporized material increases appreciably from the 
laminar value. This fraction of vaporized material is a 
function of both viscosity and type of heating (i.e., laminar 
or turbulent). For a material such as Pyrex, which has 
ablation properties similar to those of quartz except for an 
appreciably lower viscosity, the fraction of vaporized mate- 
rial is very small even for the case of turbulent heating. For 
a material of high viscosity, e.g., quartz, the fraction of 
vaporized material depends on whether the boundary layer 
is laminar or turbulent. This is illustrated in Figs. 9, which 
compare the effective energy of ablation distribution for both 


Fig. 5 Integrated effective energy of ablation around body for 
quartz. Typical 5000-nautical mile trajectory. Ry = 1 ft, 
: Reyi/Ry = 1.5, ae = 30 deg 


@ #1, A/8>>1 


«21/2, 4/8 >> 1 


*1/2,A/8 <<! 


«20 


Fig. 6 Effect of radiation on the effective energy of ablation for 
quartz. Vr = 20,000 fps, 60,000-ft altitude. Ry = % ft, 
Reyt/Ry = 3.0, a = 20 deg 


laminar and turbulent flow for the case of zero radiation. 

Figs. 9a and 9b show that the vaporization for turbulent 
flow is 50 per cent higher than for laminar flow. Fig. %c 
compares the effective energy of ablation. Figs. 8 and 9 
indicate that the effect of mass injection on Arr is small for 
the ablation of quartz, and that it is about the same for both 
laminar and turbulent flow. In this latter case, the increased 
vaporization compensates for the weaker effect of mass in- 
jection on turbulent heating rates, when they are compared 
with the laminar case at the same value of the blowing 
parameter. 

For the flight conditions indicated in Fig. 6, the tangential 
velocity at the gas liquid interface reaches a maximum value 
of about 10 fps near the sonic point. The Reynolds number 
Uw5/v» also reaches a maximum value of 3.0 at this station 
(for gases, u5/v is of the order of 10°), and the liquid layer 
thickness 6 reaches a maximum value of 5 X 10~ ft on thie 
cylindrical section of the body. 
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Fig. 7 Effects of radiation on surface temperature for quartz. __ 
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Fig. 8 Mechanisms of heat absorption. 


e = Vr = 20,000 


fps, 60,000-ft altitude. Ry = 4 ft, Rcoyt/Ry = 3.0, a. = 20 deg 


Conclusions 


The Bethe-Adams theory of ablation (10) at the stagnation 
point of a blunt body of revolution has been extended to in- 
clude the effect of turbulent heat transfer and radiation 
effects. The results for an ICBM nose cone indicate the im- 
portant result that the fraction of vaporized material for a 
highly viscous liquid layer increases appreciably for turbulent 
boundary layers—a very fortunate result for practical appli- 
cations. Also, the turbulent mass injection effect on the 
magnitude of the energy of ablation is at most of the order 
of 20 per cent, which is again a fortunate result because of the 
lack of complete theoretical understanding and limited amount 
of experimental data available on this mass injection process. 


Appendix - 
Physical Properties of Quartz (1) vet. 
Specific heat c,, Btu/Ib R 0.25 


Thermal conductivity k, (Btu/ft see R) 
3.2 x 10-4, 7 > 1800 R 
Density, lb/ft? 140 
Viscosity u, lb/ft sec 0.0672 exp[124,000/7 (deg R) — 20] 
Heat of vaporization h,, Btu/lb 5500 
Vapor pressure law 


(*;)[o4 + 0.6 (2)]" =1 
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The pressure through the gas boundary layer is 
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Fig. 9a Mechanisms of heat absorption for laminar flow. 
e= 0, Vr = 23,000 fps, 100,000-ft altitude. Ry = Wft, Royt/Ry 
= 3.0, ae = 20 deg 
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Fig. 9b Mechanisms of heat absorption for turbulent flow. 
€ = 0, Vr = 23,000 fps, 100,000-ft altitude. Ry = 14ft, Rcoyt/Ry 
= 3.0, ac. = 20 deg 
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Fig. 9c Effect of laminar and turbulent heating rates on the 
effective energy of ablation. « = 0, Vr = 23,000 fps, 100,000-ft 


altitude. Ry = ft, Royt/Ry = 3.0, ae = 20 deg 
where 


po* = exp[—104,000/7(deg R) + 18.48] 


The effect of aerodynamic pressure on the vapor pressure 
law, as given by this equation is obtained as described in the 
following. 


Vapor Pressure Reactions for Silica 
S:Oscondensed 2 S:Oogas 


S:Ocgas 2 + 30. 


Dair + De = D [A-3] 
Sle 


it 
| 4 
| 
| 
f 
= 
1.0 
uf! 
BTU ‘ 
Ib 
\- 
= 
[A-2 
: 
/ 


= 
po,!? /ps,o, 


= 0.2 Pair aa po," 


where po,’ represents O» particles created by reaction [A-2]. 
The conservation of species yields— 


pa;o, + ps;o 1 
or, since ps;0/ps;o,< 1 
me 
which with Equation [A-4] yields q 
= 3po,’ = $ psio [A-8] 
From the reactions ad and A-2] 
[A-9] 


= Ky [4-10] 


ee that po, and not poe’ is used in [A-10]. This means 

that the presence of O» from air affects the equilibrium con- 

centration of O. and hence production of S;O (Eq. [A-6]). 
Introducing Equations [A-8, A-3, A-4 and A-5] into [A-9 


and A-10] 
p[0.2p + 0.133p,.]'/? = a [A-11] 


where 


(1/2)DE(O.) — L — DE(S,02) 
kT 
For the case when there is no air, p» = p.* and po, = po,’. 
Hence [A-11] yields 
= (2/3V3)(p.*)*/# 


which when combined with [A-11] yields 
Do p 1/3 
+ 
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Nomenclature 


Cp = heat capacity per unit mass 


Cx» = concentration of ablated material in boundary layer 

fn = u/u, 

g = H/H, 

h = enthalpy (energy per unit mass) 

hess = effective energy absorbed per unit mass of ablated 
material 

l = ate <= - 

M = molecular weight ratio (air to vapor) 

n = viscosity parameter 

p = pressure 

q = heat transfer rate per unit area 

qr = 

q heat transfer rate to hot nonablating surface 

r = body dimension measured normal to axis #, 

R = body nose radius 


= time, sec 
= temperature 
= stream velocity 

= velocity tangent to body surface 
= velocity normal to body surface 
= vaporization velocity 

= ablation velocity 

= coordinate tangent to body surface 
coordinate normal to body surface — 
= thickness of liquid layer 
= viscosity 

= density 

= Prandtl number (cpu/k), Stefan-Boltzmann constant 
= shear stress 

= shear stress on nonablating surface 

= ratio of heat transfer with and without mass injection 
= radiation mean free path 

= emissivity 
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stagnation point 


Subscripts 


= gas-liquid interface 
= liquid property 
gas property at stagnation point 
property of glass vapor 

= external edge of boundary layer 
= differentiation with respect to z 
= differentiation with h respect to y 
cYL = cylinder 
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It has now been demonstrated that ablation techniques are capable of providing heat protection 
for ICBM vehicles during hypersonic re-entry into Earth’s atmosphere. In this paper, the applica- 
tion of ablation to satellite re-entry is considered. The four types of materials suitable for ablation 
are discussed: A material which gasifies but does not liquefy; a material which reacts at its surface 
with the constituents of dissociated air; a material which melts and vaporizes; a material which 
partially gasifies, the remaining material being removed by surface combustion. The result of a 
theoretical analysis of ablation for a typical satellite re-entry trajectory is presented together with 
some experimental data for materials of the four types, namely Teflon, graphite, quartz and phe- 
nolic-nylon. The satellite re-entry environment was simulated in an electric are driven tunnel. 
Comparison is also made between results in low heat transfer and high heat transfer test environ- 
ments to illustrate differences in behavior between steady-state and nonsteady-state ablation. 
On the basis of the analysis and experiments presented, it appears that little ablation of the mate- 
rials considered will occur during the ballistic re-entry of a blunt nonlifting vehicle from a satellite 
orbit. The selection of materials for heat protection can, therefore, be strongly influenced by other 
criteria. 


tural strength resulting from heating beyond the ablation 


KE VERAL techniques for heat protection of high speed re- 
zone, and the effect of local irregularities during ablation, are 


entry ballistic vehicles have been successfully demon- 
strated. Some questions remain as to the suitability of these confined to the surface. 
techniques and associated materials to the less severe but 3 Ablation materials should have a high resistance to 
longer re-entry phase associated with nonlifting satellites thermal and mechanical shock and be easy to fabricate in 


(Fig. 1). 

A typical re-entry satellite trajectory. utilized herein is 
shown in Fig. 2. With this trajectory and a blunt configura- 
tion (Rg = 3.33 ft), the boundary layer will be laminar over 
a large portion of the nose. A number of correlation formulas 
for the prediction of heat transfer exist (1-4).2_ The typical 
heating cycle used is shown in Fig. 3 and is based on an exact 
solution of the boundary layer equations by Scala (4). Three 
basie types of heat alleviation schemes are the heat sink 
(5,6), foreed mass transfer systems (4) (e.g., transpiration 
cooling devices) and self-regulating mass transfer systems 
(e.g., ablation). 


large sizes, 
Materials of interest may be grouped as follows: 

1 Plastics which depolymerize to a gas but do not liquefy 
(e.g., Teflon). 

2 Materials which sublime and react with the constit- 
uents of dissociated air (e.g., graphite). 

3 Materials which first melt and then vaporize (e.g., glass). 

4 Composite materials, such as reinforced plastics which 
pyrolyze and char (e.g., phenolic-nylon). 

The stagnation point ablation characteristics of these mate- 
rials are analyzed in the following sections. 


Ablation (7-13) has proven to be an attractive method for 
the thermal protection of high speed re-entry vehicles with Oo 
their associated sharp and short heat pulse. The application 100.0 I aera 
of this technique to the longer heat pulses and less severe = 


heating rates associated with satellites is considered herein. 


Ablation Materials 


Some generally desirable characteristics of ablation mate- # 


o GLIDE 
rials are as follows: g BALLISTIC 4 

| In general, gasification during ablation is desirable. 
The large amount of gas generated thickens the boundary Q,/Q, .20 
layer and reduces the rate of heat transfer. Gasification , 
products of low molecular weight enhance this effect because «19 
of their larger heat capacity and larger diffusion coefficients 
(14). GLIDE 

2 Ablation materials should have good thermal insulation 
characteristics so that the ablation process, the loss of struc- + 04/0, 05 Q;/Q_%1.0 
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Fig. 1 Typical re-entry heating cycles 
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Pure Gasification—Tefion 


The ablation of Teflon appears to occur by a first-order 
reaction in which the polymer depolymerizes into a monomer 
(15). The monomer has a very high vapor pressure, so that 
under most conditions it will flash directly to the vapor. 
During this process, about 750 Btu/lb are absorbed. In 
the presence of high temperature air, the | monomer may 
burn to form COF; or CO + F». ps 
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Fig. 2 Typical re-entry satellite trajectory 
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: ah Fig. 3 Aerodynamic heating of a typical re-entry satellite 
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Fig. 4 Variation of surface temperature with ablation rate for 
steady-state ablation of Tefion 


steady-state ablation 
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The first published theoretical analysis (16) considered 
hypersonic Couette flow over a Teflon wall assumed to be at 
uniform temperature. It has been found experimentally, 
however, that during the ablation process the material is still 
at the initial temperature a few millimeters below the surface 
and that a steep temperature gradient exists near the sur- 
face. The process is quite complex, and is described by the 
generalized heat conduction equation 


(v-VT + oT/ot) = Q, + V-KVT 1) 
where 
Q, = heat absorbed per unit volume by the action of (e- 
polymerization 
v = the apparent motion of the Teflon with respect ‘o 


a coordinate system which is at rest with respet 
to the ablating surface at the stagnation point 


Equation [1] can be simplified to obtain an approximate 
qualitative relationship between the ablation rate x and tie 
surface temperature of the Teflon 7,. Over a large porticn 
of any re-entry trajectory, the Teflon appears to be ablati:.g 
in quasi-steady state, so that O7'/ot = 0. In addition, tein- 
perature gradients parallel to the surface are small in coii- 
parison to those normal to the surface. Since the temper:- 
ture gradients are largest at the surface and Q, is an Arrhv- 
nius function of temperature, most of the heat absorption hy 
depolymerization occurs near the surface; hence, Q, can be 
assigned to the surface itself. The first-order reaction rate 
then yields 


E = 81,400 + 900 cal/mole 
B = (9.4 X 
Combination of Equations [1 and 2] a integration using 


constant thermal properties yield the following relation be- 
tween T,, 


= — To/Tw) [3] 


The plot in Fig. 4 indicates that the surface temperature is 
relatively insensitive to the ablation rate x 
This insensitivity makes the Teflon ablation problem simple 


to evaluate. A laminar boundary layer analysis (17,18), in 
_ which all species are assumed to have the same physical 
properties and the pressure gradient is neglected, yields the 
following expression for the energy transfer in the presence 
_ of mass transfer of a combustible gas 


—(K OT /0x) pw = ~ + AH Co) [4] 


—(K 0T/02) pw = 


- Since the wall temperature is insensitive to the ablation rate, 


H,, and AH, may be taken at some average value, and Equa- 
tion [6] yields the ablation rate directly. 

For mass transfer into a laminar boundary layer, Stewart 
(14) has correlated the results of theoretical analyses to ob- 
tain the following expression for the Stanton number 


10" a 0 072 (* = | 7 
X , ABLATION RATE,FPS Pelle! \ 4 


:/Cp,i,) is effectively 29/Mc,r,. If the heat transfe: . 
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given by Equations [6 and 7], is equated to the product of the 
enthalpy increase of the plastic and the steady-state mass 
transfer rate, the following expression for steady-state value 
of the heat of ablation is obtained 


AH, 
] + w)air 
With 7, taken at an average of 1340 F, JT) = 100 F, 
C,, = 0.30, AH, = 738 Btu/Ib, AH, is equal to 1110 Btu/Ib. 
The heat of ablation of Teflon is shown in Fig. 5 for which 
a value of AHy = 9380 Btu/lb was used. If combustion is 
ignored, the theoretical prediction is too large by 20 per cent 
at the lower enthalpies. At high enthalpy levels (i.e., signifi- 
cant portion of trajectory), the theory of (4) may be used 
for the energy transfer to the surface 


Q = Qn-o — 0.48 (H. — [9] 


+ 0.439(H, — Hu)air 


The transient heat conduction equation has been solved 

using Equation [9] and the method of (19) for the re-entry 

trajectory heating rate shown in Fig. 3. The resulting 

ablation rate is shown in Fig. 6. Some temperature profiles 

are shown in Fig. 7. Total ablation is 0.2in. = 


Surface Combustion—Graphite 


Graphite is typical of the class of materials which undergo 
surface combustion (20). 

In the absence of detailed reaction rate data, a simple 
model for the complex C—O—N system is employed, and 
certain simplifying mathematical assumptions are generally 
made. 

The following assumptions were made in the analysis ‘per- 
formed by Scala (21): 

| The gas phase can be approximated by a three-com- 
ponent mixture, consisting of nitrogen, atomic oxygen and 
carbon monoxide. 

2 Gas phase reactions are assumed frozen. 

3 Nitrogen reactions do not occur at the surface. 

4 The surface combustion product is carbon monoxide. 

5 The Lewis and Prandtl numbers are constant and equal 
to |.2 and 0.75 respectively. 

It can then be shown that the ablation rate of graphite at 
the stagnation point of an axially symmetric body is given by 
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Fig. 6 Material ablation rates for re-entry satellite trajectory 
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Fig. 7 Temperature profiles in Teflon after ablation ends 


This solution applies when the slowest, and hence rate- 
determining, step is the process of diffusion and counter- 
diffusion, where the ablation rate depends very little on the 
surface temperature. 

Assuming a Newtonian pressure distribution over the 
body, Equation [10] may be further simplified to the form 


Rp = 5 X 10* VP3/p, 11] 


Substitution of typical satellite re-entry pressures and 
densities into Equation [11] yields the values of m, shown in 
Fig. 6. Thus, in the case of satellite re-entry, graphite has 
such low thermal erosion characteristics that it may be con- 
sidered primarily as a heat sink material which undergoes 
oxidation reactions with the environment. 

Considering heat transfer into the interior of the graphite, 
an energy balance reveals that 
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Since the instantaneous combustion rate m, is very small 
during satellite re-entry, the ablation process should not alter 
the heat transfer to the surface appreciably. Further, since 
How actually can be a negative quantity (it contains a rather 
large negative heat of formation of the combustion products), 
it has been found that for the most severe flight conditions, 
the chemical heat release actually exceeds the heat blocking 
term owing to boundary layer thickening, so that —(K X 
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Fig. 8 Front and back wall temperatures of graphite for re-entry 
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Fig. 9 Variation of Q* with stagnation pressure and stagnation 
for fused silica 
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Fig. 10 Surface temperature history for fused silica 


d07'/dx).0 is approximately 20 per cent larger than Qaero at 
worst. 

Thus, the design of a thermal shield utilizing graphite is 
limited not by its ablation rate, but rather by the heat pene- 
tration aspects and, consequently, mechanical limitations, 
such as thermal shock and shield fastening problems. The 
front and back wall temperatures of a layer of graphite 0.10 
in. thick, with a surface emissivity of 0.8, are shown in Fig. 8, 
for the case of variable thermal conductivity. 


Melting and Vaporization—Glasses 


There are few practical materials which sublime directl\ or 
which have sufficiently high vapor pressures at their melting 
point to give flash vaporization. To obtain appreciable 
vaporization, with its associated heat absorption and bound- 
ary layer thickening, materials are required whose mel is 
highly viscous, so that the liquid can be heated to n«ar- 
boiling temperatures before aerodynamic forces remove it. 
Such materials are by their very nature glasses, e.g., liquids 
which are viscous enough so that they may be supercocled 
without having their molecules falling into the ordered airay 
of a crystal. Data on glasses indicate that fused silica jias 
the highest viscosity, and under certain conditions ‘his 
material has given evidence of superior ablation characier- 
istics (7,11). 

In addition, high silica content glasses are quite resistan' to 
thermal shock. They are, however, transparent to infrared 
radiation, and consequently often have additives dispersed 
throughout to absorb or scatter both incoming radiation «nd 
self-radiation from the heated surface. 

When vaporized, silica forms both the gaseous dioxide and 
monoxide, in the following highly endothermic reactions 


With the formation of a liquid phase, the effect of body 
forces owing to deceleration of the vehicle during re-entry 
must be considered. The body forces can slow the liquid 
sufficiently so that it detaches from the surface somewhere 
downstream, leading to a decrease in Q*. 

Fanucci and Lew (22) analyzed the flow of a liquid layer 
with vaporization under the influence of a gaseous hypersonic 
boundary layer (laminar or turbulent, depending on the 
region of the body concerned), taking proper account of the 
history of the flow from the stagnation point. Some results 
of the characteristics of the liquid phase at the stagnation 
point of the satellite vehicle have been matched with an anal- 
ysis by Scala (23) for the dissociated hypersonic boundary 
layer with mass addition to obtain the rate of ablation of the 
surface, the fraction vaporized [ and the effective heat of 
ablation Q*. 

The heats of ablation for steady-state conditions are shown 
in Fig. 9. An increase in stagnation pressure causes the sur- 
face temperature to rise, resulting in a decrease of the vis- 
cosity of the liquid layer and a decrease in Q*. The lessened 
viscosity causes more material to flow downstream under 
the influence of aerodynamic forces, resulting in a smaller 
gasification ratio. On the other hand, Q* increases in a 
nearly linear fashion with increasing stagnation enthalpy 
because the higher heat transfer rates at high stagnation 
enthalpies force more of the liquid to vaporize. 

The variation of Q* along the trajectory of the re-entry 
satellite is indicated in Fig. 9. The associated surface tem- 
peratures are shown in Fig. 10. There is little surface tem- 
perature variation for steady flight at each altitude and 
velocity of the re-entry satellite trajectory. It is clear that 
during the actual trajectory the initial temperature of the 
silica is low, and no ablation will occur until the silica is heated 
above its softening point. This occurs at about 160 sec 
(see Fig. 3) when the ablation process begins. 
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A study of the effects of deceleration, and hence of body 
forces, indicates that Q* is increased and the net mass loss 
over a trajectory is decreased. This effect is shown in Fig. 
6 and is not negligible even though the rate of deceleration is 
low. Surface reradiation will decrease the net heat transfer 
to the surface. In Fig. 6 it also may be seen that a very 
modest emissivity of 0.1 decreases the ablation rate by about 
20 per cent. If the surface emissivity is unity, then the 
equilibrium temperatures will be 3470 R for the case of no 
heat conduction into the interior. Under these conditions, 
the silica probably will not ablate but will act merely as a 
reradiating surface. If the thickness is sufficiently large, 
then the silica will act as a heat sink, with a maximum sur- 
face temperature of 2800 F. A heat sink calculation yields 
a maximum value of the net heat transfer of 28 Btu/ft?-sec 
and a heat pulse of about 5000 Btu/ft?. The thickness of 
silic. required to yield a final back face temperature rise of 
200 F is about 7 in., and is inversely proportional to the final 
back face temperature rise; a thickness of 1 in. will yield 
a back face temperature rise of 1400 F. On this basis, it ap- 
pears more desirable to have the emissivity of the silica less 
tha unity so that a surface temperature sufficiently high to 
cause ablation is reached, resulting in a decrease of the net 
heat content of the remaining silica. Thus quartz, when 
utilized as a heat sink, is relatively inefficient. 

Nonequilibrium effects have been considered by Scala and 
Vidale (24,25). It was found that rate controlled vapori- 
zation tends to increase the surface temperature resulting in a 
decrease of T and Q* above 200,000 ft. For the vehicle and 
the trajectory shown in Fig. 2 the effect is not significant, 
since ablation probably will not begin at altitudes above 
200,000 ft. 


Reinforced Plastics 

Reinforced plastics include some of the most attractive 
materials for use in satellite re-entry. Certain of the phe- 
nolic, silicone and melamine resins are suitable plastics 
(7); 91LD resin (Cincinnati Testing Laboratories), for ex- 
ample, has shown good characteristics experimentally. 
When heated, this resin pyrolyzes into a gas and a carbonace- 
ous structure which shrinks in dimension. Preliminary 
experiments indicate that the gas is composed of hydrogen, 
methane and carbon dioxide. Combustion may occur be- 
tween methane and oxygen to form CO:, H.O, OH and CO. 
The latter three then may react with and chemically erode 
the carbonaceous structure to form CO and CO». Suitable 
reinforcements include “E”’ glass, Refrasil (a form of silica: 
H. T. Thompson Co.) and nylon (7,11). 

Phenomenological behavior in this class varies signifi- 
cantly with time and magnitude of heating. For the appli- 
cation under consideration, the nylon reinforced plastic 
pyrolyzes to form a moderately thick, low density, charred 
layer which, in general, has shrunk dimensionally and which 
may have been chemically eroded at the surface. The surface 
temperature of the char layer may become sufficiently hot 
to reradiate most of the aerodynamic heat. The heat 
damage thus propagated into the material is caused by the 
aerodynamic heat transfer and by the heat released during 
— or chemical reaction within the char and adjacent 

yer. 

For glass and Refrasil reinforcement, char depth was 
found experimentally to be about twice that of the nylon 
reinforcement. When inorganic fibers are used as reinforce- 
ment, the plastic tends to delaminate and grow in dimension 
(13). The local rate at which unpyrolyzed plastic is pyro- 
lyzed is given by 


Opr 
The total rate of pyrolysis is given by dar 27-67 
f, PrpBe [15] 


1960 


The local rate of heat absorption is given by | 


Qe = pop AH e~2/RT (16) 


Thus the problem of predicting the rate of pyrolysis of the 
plastic is quite complex. Steady-state ablation does not ap- 
pear probable for re-entry satellites protected by reinforced 
plastics. Therefore, redefining Q* by Q* = Qcm =0)/my, where 
My is the rate at which mass is pyrolyzed but not necessarily 
removed, Scala has obtained the expression 


Q* = + ~H» [17] 


To evaluate Q*, the chemistry of the interaction between 
the material and dissociated air must be considered. 

If pyrolysis occurs in a vacuum, then I'm, represents the 
rate of gasification, and 1 + (1 — I')m, represents the rate of 
formation of char, where m, is the rate of pyrolysis. Dur- 
ing hypersonic re-entry, however, oxygen diffuses toward 
the surface and oxidation reactions occur. The rate of gasi- 
fication is then given by (Tf + =)m,, and the instantaneous 
rate of char formation is (1 — I — =)m,. Char is con- 
tinuously formed at an interface between the virgin plastic 
and the pyrolyzed layer, and, if sufficient oxygen diffuses to 
the surface, is consumed at the outer surface. Thus the rate 
of growth of the char layer depends on: 

1 The mass diffusion flux of oxygen from the main 
stream toward the surface. 

2 The quantity and reactivity of the combustible gases 
produced during pyrolysis which flow into the boundary 
layer; this, in turn, is controlled by the rate at which heat 
can be conducted to the virgin plastic. 

For the case of ablation of phenolic nylon, an exact solu- 
tion for air injection into the hypersonic dissociated stagna- 
tion point laminar boundary layer yields the following ex- 
pression for the heat blocked due to mass transfer (23) 


Am 


4 


(2) [240 + 0.48(H. 


This is assumed to be valid for predicting the effects of gas 
injection containing species other than oxygen and nitrogen. 
Since the mean molecular weight of the gas at the surface 
appears in the denominator, the injection of lightweight gases, 
such as hydrogen, enhances the blocking term. 

If nylon is utilized as a resin, it melts, vaporizes and reacts 
without leaving a char residue (11), whereas phenolic-resin 
decomposes to form a char residue. Nylon and phenolic- 
resin combined in equal amounts by weight may be repre- 
sented by the empirical formula C,3H,O.N. The behavior 
of the composite material is such that lightweight gases are 
released during pyrolysis and a char layer is formed which 
enhances reradiation from the surface. 

A theoretical analysis of the ablation of phenolic-nylon, 
utilizing the boundary layer equations, is complicated by the 
large number of different chemical species which are the 
original products of pyrolysis and by the subsequent com- 
bustion reactions between these species and the components 
of dissociated air. Experimental data (11) indicate that 
under vacuum conditions T & 0.7. Hence, 0.7 < (f + 2) 

The dependence of (AQ/Am),, the heat blocking action 
term, on the mean molecular weight M,. has been noted in 
the foregoing. If, for example, the final products of gaseous 
dissociation and combustion are primarily hydrogen, with 
small amounts of water and carbon monoxide, the heat 
blocking term can easily be a factor of 2 or 3 larger 
than the corresponding term for transpiration cooling by air 
injection. Another effect contributing to Q* is the magni- 
tude of the gaseous enthalpy at the surface H,. Combining 
Equations [17 and 18] indicates that Q* tends to increase 
approximately as the additive term (3)H.. The more — 
hydrogen at the surface, the larger the term H, will be. 
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Finally, it is observed that complex hydrocarbons, such as 
phenolic-nylon, have a negligible heat of formation per unit 
mass (Hyp ~ 0). 

Evaluating Equation [17] in the light of the foregoing dis- 
cussion, Q* may be approximated by (26) 


“ail 
Q* = 1500+ (1/2)H. [19] 


In high speed re-entry mass loss due to ablation (related 

to Q*) is generally the primary consideration; in satellite 
re-entry the temperature reached at the inner wall surface is 
perhaps of equal interest. Since conductivities for all of the 
materials considered here* (except graphite) are of the same 
order of magnitude as the conductivities of the best practical 
insulators, the self-insulation thickness is a good figure of 
merit. 

In view of theoretical complications introduced by long- 
time heating (e.g., charring and the prevalence of non- 
equilibrium phenomena), experimental data are of par- 
ticular interest. 

Experiments were performed in a Mach 4.8 arc heated 
supersonic wind tunnel at a stagnation enthalpy of 150 RT», 
calorimeter heat transfer rate of 65 Btu/ft? sec and a stag- 
nation pressure of p. = 0.026 atm.‘ Tests were run for 
1 and 2 min; the total heat pulse for a 2-min run is comparable 
to that shown in Fig. 3. 

Thermocouples were placed at two positions in the models 
(Fig. 11) so that the weights of material between the thermo- 
couple and initial surface were 1.15 lb/ft? (location A) and 
2.30 lb/ft? (location B). The thermocouples were at loca- 
tions A and B for the 1-min runs and at B for the 2-min run. 

Data obtained at the end of 1 min (thermocouple at A) and 
2 min (thermocouple at B) are shown in Fig. 12. The 
temperature rise at each station at the end of each run is re- 
ferred to the temperature rise of graphite. It is seen that 
for both heating intervals, the minimum temperature rise 
occurs with phenolic-nylon. 

In view of the fact that the time from the end of significant 
heating until impact may be several times as large as the time 
of heating, temperatures have been recorded after as well 
as during heating. These data are shown in Fig. 13 for sta- 


3 Graphite, types ATJ and ATL; phenolic-nylon, 40 per cent 
resin by weight; phenolic-Refrasil, 38-42 per cent resin by 
weight; phenolic- glass, 40-50 per cent resin by weight. 

he otal carbon contamination approximately 5 per cent by 
weight. 


820 


THERMOCOUPLE 


(0 STA 1- HEAT ON {MIN 
STA 2- HEAT ON 2MIN. 


2 
az oun = 
- wd <a 
Cia * 
tall 
g}— 
lac: 
|S 


WwW 
Bes REFRASIL 
a6 
= 
WwW 
PHENOLIC - 
SS 


GLA 


PHENOLIC - NYLON 


! 2 3 4 5 6 7 8 9 10 
TIME — (MINUTES) 


' 
HEAT OFF 


Fig. 13 Temperature histories of material specimens after heat 
cutoff 


tion B for a 1-min initial heating load. It should be noted 
that heat loss in the tunnel under these conditions is pri- 
marily by radiation, so that the convective process and 
combustion occurring during actual re-entry are not simu- 
lated. Although both of these processes probably act in 
opposite directions, only qualitative conclusions are drawn 
here. It is seen that, with the exception of graphite, the 
maximum internal temperature occurs after the heating 
cycle has been terminated. 

The change in position of the temperature rise of Teflon in 
relation to phenolic-nylon at the end of the 1-min run at 
stations A and B is of interest (Figs. 12 and 13). The higher 
temperature of the Teflon at station A probably is due to the 
greater ablation rate of Teflon in these tests, which brings 
the surface of the Teflon quite close to station A. The lower 
temperature of the Teflon at station B indicates that the 
temperature gradient between station A and station B is 
greater than that in the case of the phenolic-nylon. Thus, 
as would be expected, less energy is stored in a material which 
has a higher ablation rate during the heat pulse. The high 
surface temperature material will, however, reradiate more 
of its stored heat after the pulse, so that the penetration of the 
shield material by heat is reduced. From the data shown 
in Fig. 13, it appears that this latter effect is responsible for 
the lower, late-time temperatures in phenolic-nylon. It is 
possible that the trends shown for Teflon and phenolic-ny!on 
could change for different heating rates, total heating time, 
shield weight and thickness. 
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Fig. 14 Model cross sections—heat soak study: A, phenolic- 
Refrasil; B, phenolic-nylon; C, phenolic-glass. p. = 0.026 
atm; H./RT, = 150; Q,.. = 65 Btu/ft?-sec 


Table 1 Measured heats of ablation 
* * 
(theoretical), (experimental), 
Material » Btu/lb Btu/lb 

Teflon ; 2700 3120 
phenolie-nylon «4180 4480 
phenolic-Refrasil 7530 
phenolic-glass 9600 


Corrected for surface temperature. 


Test Conditions 


Quo = 65 Btu/ft®sec 
H./RT, = 150 
Pe = 0.026 atm 


The measured heats of ablation are shown in Table 1. 
Although the last two materials exhibit large values, it is 
not evident that they would be the preferable because of the 
heat conduction problem, particularly in the case of graphite. 

The measured heat of ablation of Teflon is somewhat 
higher than the value obtained during high heat flux condi- 
tions (Fig. 5), probably since heat is conducted into the 
material faster than it ablates, resulting in 


mAH, < (K 07/02) pw 


Thus, the measured heat of ablation in this test should be 
larger than the steady-state value in which 


mAH, = (K OT /02) pw 


Nonsteady conduction could also account qualitatively for 
the value of Q* for phenolic-nylon, which is larger than the 
theoretical prediction. 

Sections of the models tested are shown in Fig. 14. The 
formation of a thin layer of solidified glass may be observed 
on the surface of the phenolic-glass model during the latter 
portion of the 2-min heating cycle; thus basic phenomeno- 
logical differences in the ablation process with time are to be 
expected. 

The significant variation in behavior with heating rate for 
some nonhomogeneous materials is shown in Figs. 15 for 
melamine glass. Total heat pulses were essentially identical; 
maximum heating rates differed by a factor of 14. Major 
differences in characteristics of the char layer are clearly 
evident. In the case of the low value of maximum heat 
transfer rate, an extensive fibrous char structure of pyrolyzed 
plastic exists, covered on-the surface by molten reinforcing 
material. In the case of the high value of the maximum 
heating rate, the higher ablation rate clearly limited the 
depth of the remaining char layer, resulting in a higher mass 
loss, but a relatively undisturbed remaining structure. 
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TEST CONDITIONS 


175 = 175 
Q= 60 BTU/FT?SEC; t=I20SEC Q=830 BTU/T SEC, t=10 SEC 


2 
7200 BTU/FT = 8300 BTU/FT? 


TEST CONDITIONS 


175 175 
Q= 830 BTU/FT°SEC; t=I0 SEC Q=60 BTU/FT*SEC; t= 120 SEC 
= 8300 BTU/FT® = 7200 BTU/FT? 


TOTAL TOTAL 


Figs. 15 Comparison of material behavior at high and low abla- 
tion rates 


-, 


Conclusions 


On the basis of the analyses and experiments discussed, it 
appears that very little ablation of the materials considered 
will occur during the ballistic re-entry of a blunt nonlifting 
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vehicle from a satellite orbit. The largest amount of abla- 
tion will occur with Teflon, but will amount to less than 3 Ib 
per ft? for the re-entry heating cycle considered. 

The selection of the materials for heat protection can there- 
fore be based on other criteria. Foremost of these is the 
weight of material required for self-insulation, since it appears 
likely that the ablation material also can serve effectively as 
the thermal insulator. On this basis, both Teflon and 
phenolic-nylon qualitatively appear to be the most suitable 
of the materials investigated. If we also consider structural 
integrity, the phenolic-nylon has advantages over Teflon, 
since the recrystallization of Teflon and the concurrent ex- 
pansion present serious structural problems. 

On the basis of conservative estimates, an initial Teflon 
thickness of 1 in. would be more than adequate to maintain 
the inner surface temperature rise to a value below 200 F. 
The same weight of phenolic-nylon would have equivalent 
performance. The best available heat sink material would 
require not only a larger weight to absorb the total heat 
transferred to the body, but also a considerable weight of 
insulation material to protect the payload from overheating. 
It is probable, therefore, that an ablation heat shield struc- 
ture will be approximately half the weight of that required 
for the best heat sink material. Presumably this ratio can 
be further improved by utilization of suitable lightweight 
insulation in place of the self-insulation thickness of the 
ablation material. 


M = molecular weight 
= mass transfer rate = ppé for steady-state ablation 
Pp = pressure 
a ie = heat transfer rate to surface 
Q* = heat of ablation = Q (m = 0)/my 
Qim=o) = heat transfer rate in the absence of mass transfer 
p = energy absorbed per unit volume due to pyrolysis 
(AQ/Am)wo = change in energy transfer per unit change charge 
in interphase mass transfer for no combustion 
= universal gas constant 
oll = nose radius of the re-entry satellite 
= time 
temperature 
tangential velocity parallel to surface 
velocity vector 
wae coordinate normal to surface, direction into the 


material 


= ablation rate 


= thermal diffusivity 
= vaporization coefficient _ 
= stagnation point velocity gradient 

= fraction of material ablated which is gaseous 
€ = total surface emissivity 

viscosity 


= fraction of virgin plastic that enters gas phase 


Nomenclature 
B = frequency factor for depolymerization : 
Cx(m) = Stanton number in the absence of combustion 
Cy = concentration of oxygen ; 
C> = specific heat at constant pressure 
E = energy of activation of depolymerization 
H = gaseous stagnation enthalpy, including chemical 
enthalpy 
AH, = heat of combustion per unit mass of oxygen 
AH, = enthalpy change, initial to surface fempersiece of 
w 
plastic = Hpw(g) — Hpo(s) = Cp,(s)dtAH 

AH, = heat of depolymerization per unit volume for 

2 plastic; latent heat of vaporization y 
K ——_— = thermal conductivity 


p = 
Subscripts 
aero = aerodynamic 
e = outer edge of gaseous boundary layer, stagnation 
conditions 
- 
= gas. 
0 = initial 
Pp = plastic 
t = total 
vp = virgin plastic _ 
w = wall 
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nals of the Ablation of 
Plastic Heat Shields That 


Research and Adv anced Development Div., — ' 


3 The ablation of plastic heat shields that form a charred surface layer under extreme aerodynamic 
5 heating is analyzed. The model chosen for this type of ablation consists of a porous outer layer 

— ‘through which there is a counterflow of gases from the thermal decomposition of the plastic. The 
ss @harred layer is backed by a semi-infinite layer of the undecomposed solid. The transpiring gases 
-. { are considered to be in thermal equilibrium with the skeleton of the porous material. The heat 
aft transfer equations are set up, and the steady-state solution obtained. The formulas for the steady- 
_—s- state temperature distributions, charred layer thickness and heat of ablation in terms of model 
-_- parameters and variables of the thermal environment are derived. An approximate formula is 
also derived for the induction time of the nonsteady phase of ablation. 


Description of the Model tions of the two surfaces, which are given by 
N DER conditions of extreme aerodynamic heating, heat 8, = f v,(t) dt 
shielding may be accomplished by consumption of the ° 
shielding material (1). This is known as ablation. During 


where v,(¢) and v,(t) are the velocities of the front surfaces of 
the solid and the layer, respectively. 
We will assume that the charring process is such that 


the ablation of many materials containing organic matter, 
such as plastics or resins, a porous charred layer is formed on 
the surface (2,3). This layer may attain a thickness of 
several millimeters in the quasi-steady state. It complicates 
the heat transfer problem not only because it introduces a 
second phase between the original solid and the boundary 
layer, but also because it contains a counterflow of gases from aero ae jee FEN Mere: 
the decomposing solid. 1 

In addition to the usual outside heat transfer surface next 
to the aerodynamic boundary layer, we have in this case 
another boundary to consider between the charred layer and 
the original solid where a sharp change in properties occurs 
as well as an absorption of heat AH- by chemical decomposi- 
tion. We can therefore expect a discontinuity in the tem- 
perature gradient at that second interface. During ablation, 
both interfaces will be moving. In the transient phase, the 
two interface velocities will be unequal. 

The geometry of the model and the temperature distribu- 
tions are illustrated schematically in Fig. 1. The frame of 
reference 7’, X is considered fixed relative to the material of 
the original uncharred plastic. Another frame of reference 
T”,, X", is specified to move with the front surface of the Flow 
charred layer at = A third frame of reference 
X’, is specified to move with the back surface of the charred = 
layer (i.e., the front surface of the original solid) at z = s,. a 
We therefore have the relation cae 


L = 8 — 8 (1] 
between the thickness of the charred layer L and the loca- - bi 
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° Senior Scientist, Theoretical Physics Group. Fig. 1 Schematic illustration of the model for ablation of plas-— 
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voids are formed in the plastic in such a way that there is Cc = heat capacity 
no overall shrinkage and that the material points in the two K_ = thermal conductivity 
solid phase structures do not move relative to each other. q(x) = heat flux in the positive direction at the station .c 


We can therefore use a common, fixed frame of reference for 


the two solid phases. One of the consequences of this as- K, C and p are assumed to be constants. In this frame of 


reference, ablation is represented by the motion of the 


a ; boundary condition at the surface of the solid g(x) which 
mM, = m.(p. — pi)/ps a [4] must allow for heat absorbed by any phase changes tliat 
may occur. 
oe ~~ oe [5] In later discussions, it will be useful to consider the hat 
t- & We also need to make some assumptions about the structure flux at z in terms of the increase in enthalpy H in the region 
and properties of the porous charred layer and about the to the right of a fixed vertical surface at x 7 
nature of the gas flow and its interaction with the solid skele- i. t dH 
ton of the layer. These assumptions are that the size of the H = pC f. T(x)dz = / ad dt 9] 
voids which constitute the porosity of the charred layer 
is sufficiently small, compared with the charred layer thick- Since enthalpy must be conserved, the rate of increase ii H 
ness, so that the counterflowing gases are always at the same must equal the heat flux across the surface 
temperature as the solid skeleton at corresponding points 
where they are in contact, regardless of the magnitude of ofa) dH _ oC OT (x) da 10) 
the temperature gradient. Weinbaum and Wheeler (4) have dt 0 ot 


analyzed the problem without the latter assumption and 


showed that the gas temperature rises very quickly (within Heat Transfer Relations in a Frame of Reference That 


a layer about one pore diameter in thickness) after it enters ; ; : 
at the back of the porous structure. We also assume that ee 
the gas flow is unimpeded by the porous layer; this means In Fig. 2, the coordinate system T’,X’ is moving with re- 
that the flow rate and gas pressure are not functions of the spect to the stationary system 7',X (and with respect to the 
coordinate z. We will neglect heat transfer by radiation material of the body) at a speed v(t) which need not be con- 
across the voids and by back diffusion in the gas stream. stant. The relationship between coordinates in the two 
Finally, we will assume that the thermal and physical frames of reference is therefore given by the two equation 
properties of the three phases are constants, independent of : 
temperature. We will, in particular, define the properties 
of the porous layer, such as its density :, thermal conduc- a ee t sisal 
tivity and heat capacity to include the gas enclosed in o(t)dt [12] 


its pores under conditions of no gas flow. 

It is the assumption of constant properties in particular 
that make the model an idealized one. These properties are 
expected to vary with temperature and/or depth in the 
material (3). For engineering estimates of ablation be- 
havior, however, it will be justified to use the formulas so 
derived with average or equivalent-model values of the 
basic parameters. Moreover, this ideal model can be ex- OT(x) dT"(zx’) oT (zx) 
pected to demonstrate the essential characteristics of real ogee —v(t) “ 
material ablation in a general way if not in detail. 


The temperature is the same at corresponding points in 
the two coordinate systems at all times. Therefore all 
slopes and curvatures are the same at corresponding points 
in the two systems. The time derivatives of T are different 
in the two systems and are related by the usual equation for 
the convective change (6) 


[13] 


: The equation for the convective change, when applied 
the heat flux as defined by Equation [10], leads to the 
Ablation Without Charred Layer Formation relationship between the heat fluxes in the two frames of 


The heat transfer relations for a body that ablates without _ Teference : 


forming a liquid or charred layer are known (5). They will q'(z’) = q(x) — pCv(t)T (x) AY 
be reviewed here since they form part of the solution of the : bey 
present problem. The discussion will also serve as an in- on ~~ ” 
troduction to the important concepts of ablation used in the + 

T 


analysis of our model. 


Heat Transfer Relations in a Frame of Reference That Is 
Stationary With Respect to the Material of the Body 


In Fig. 2, the coordinate system 7, X is stationary with re- 
spect to the material of the body. The conventional heat 
transfer equations for a one-dimensional case are applicable 
in this coordinate system 


7 = temperature increment above initial temperature 
t = time 
p = density 


Fig. 2 Temperature distribution in heat transfer system in 
two frames of reference moving with respect to each other 
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pCv(t)T(2) is the rate of enthalpy loss across the surface at 
y’ in the T’,X’ coordinate system due to the material’s 
relative backward motion. 

In order to obtain the heat transfer equations in this co- 
ordinate system, one merely needs to substitute Equations 
(13 and 14] in Equations [6 and 8] with x = 2’. The form 
of the boundary conditions will be discussed in the next 
section. 


Ablation With Vaporization Under Aerodynamic Heating 
Conditions 


In the ease of vaporization during ablation under aero- 
dynamic heating, the net aerodynamic heat input is a bound- 
ary condition that moves with the ablating surface and can 
be related to q’(0) which is given by (7,8,9) 


q'(0) = q(aerodynamic) — q(transpiration) — 
g(radiation) — m.Hr [15] 
where 
g aerodynamic) = CyopeU — hs) [16] 
y(transpiration) = mon(h. — hs) [17] 
Mm, = pv(t) [19] 


= gas density at edge of boundary layer 


n = constant of the gas 

m = mass flow or ablation rate 

Hr = latent heat of vaporization plus the 


sensible enthalpy of solid at 7’(0) 


y(transpiration) represents the “transpiration” or “blowing”’ 
effect, which is the amount that the aerodynamic heating 
rate is reduced by the countercurrent of gases from the va- 
porization of the solid. This is a complex phenomenon that 
has been the subject of much theoretical analysis (7). Rigor- 
ous solutions of the boundary layer equations are very cum- 
bersome and require extensive material properties that are 
usually not available. Fledderman (8) has suggested Equa- 
tion [17] is an approximation of sufficient accuracy for simple 
engineering analysis. The constant 7 varies approximately 
as \f-'/3 or M-/4, where M is the molecular weight of the 
gas. 

From Equations [8 and 14], we can write the boundary 
condition at the surface for the heat transfer equations in the 
stationary frame of reference as 


(Cuop.U, — mon) (he — hs) — oe[T(s)]4 + 


pCr(t)T(s) — m.Hr [20] 


The sensible enthalpy was included in Hr to keep Equation 
[20] in a form related to Equation [14]. There will be a 
cancellation of sensible enthalpy between the last two terms 
of Equation [20]. This complication arises from the intro- 
duction of a phase change into the model. 

Because of Equations [13 and 14], we can write the dif- 
ferential equations for the temperature distribution in the 
moving frame of references as 


= —K — pCo()T"(x') 


The upstream boundary condition is 


= K }. (Cxopel m.n)(h. hs) 
lo 


+ p Co(t)T’(0) — m.Hr [23] 
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Quasi-Steady-State Relations 


For the quasi-steady-state solution, we let 07’/ot = 0 in 
Equation [21], integrate twice and evaluate the constants of 
integration from the boundary conditions 

T,(0) = surface temperature [24] 
(Subscript * designates steady-state condition.) 
T,(2) [25] 
The result is the familiar relation (5) 
T(x’) = T,(0) exp (—2"/d) [26] 
The characteristic depth d is given by 
+, 


A = K/v.apC 


From Equation [10], we therefore get 
qx(s) = pCr.xTx(s) = Cm Tx(s) [28] 


and from Equation [14] we get 
q'x(0) = 0 [29] 
The heat stored in the body in the steady state is 
Hs = pCdrT,(0) = T,(0)K/vse [30] 


from Equation [9]. 

These results indicate that in a practical case the heat 
flux into the surface of the body is initially equal to the net 
heat flux to a cold surface in either frame of reference. The 
heat flux into the body in the stationary frame of reference 
q(s,t) decreases asymptotically to cm.«7*(s) as time goes on 
and the steady state is approached. The heat flux into 
the body in the moving frame of reference q’(0,t) asympto- 
tically approaches zero as time goes on. 


Induction Time for the Nonsteady Phase 


The steady state of ablation represents a limit that is ap-— 
proached asymptotically, and the time to reach the steady 
state is therefore infinite. The characteristic time 7’ for 
this transient nonsteady phase is sometimes taken as the 
ratio of the heat stored in the body to the heat input rate (5) 


= \pCTx(s)/q(0,0) [31] 


There is another quantity, which we will call the induction 
time 7, that has considerable practical importance. It is 
the intercept, on the time axis, of the extrapolated steady- 
state position-vs.-time curve of the ablating surface. This 
is illustrated in Fig. 3. 

The two quantities 7’ and 7 are equal in the special case 
where radiation losses from the surface and hot wall correc- 
tions for the aerodynamic heating are neglected. In this 
case, the total amount of heat q-* carried away from the 
body per unit mass of ablated material is constant for all 
time. The total sensible heat H per unit area which has 
accumulated in the body at any time ¢ is then equal to the 
total heat flux q(0,0)t into the body from the time heating 
started minus the total heat pq-*s which has been carried 


away by ablation. Thus 
H = q(0,0)t — pq-*s [32] 
After sufficiently long times, s is given by 
S«(t) = v«(t — 7) [33] 
and 
H = q(0,0)t — pqc*ve(t — 7) = q(0,0)r [34] 
since 


q(0,0) = pqo*ve [35] 


€ 
3 
ae 
4 


0 Tt 


os 3 Location of ablating surface vs. time 


Thus the induction time 7 is adi the time which would be 
required to accumulate the heat stored in the steady-state 
temperature distribution of the body, if no heat were carried 
away by ablation. 


Heat Transfer in the Porous Charred Layer 

The heat transfer equation for the porous charred layer 
of the model we specified earlier can be written (in the 
stationary frame of reference) simply as (10) 


The second term on the right-hand side of the equation is 
simply the convective term for the gas flow. The heat 
capacity and conduction contributions of the gas are in- 
cluded in C; and Ki. 

The boundary conditions in the stationary frame of refer- 
ence for the layer in our composite system (i.e., the layer 
plus the semi-infinite solid in back) are 


qi(si,t) = Ox1)2,=s, [37] 


Ti(21,0) = 0 [38] 
— = qs(8s,t) + m, AH. [39] 
where 


qu(si,t) = heat flux in stationary frame of reference at 
front surface of the layer 
 qs(8s,t) = heat flux at the front surface of the solid, i.e., 
*hy heat flux into the solid body in the stationary 
; frame of reference 
AH. = heat of charring the solid to form the gas and the 
layer material 


For the boundary conditions at the back surface of the 
layer (Eq. [39]), we assume that the isothermal -decompo- 


sition reaction 


golid ——> gas + layer [40] 
T’.(0) 


t) = [42] 
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occurs; and that the layer material is generated at the 
decomposition temperature (i.e., the front surface tem- 
perature of the solid). The consequence of this is that 
T"(L) = T'(0). Therefore, the heat flux coming into the 
back surface from inside the layer is accounted for by A//, 
and further heat transfer into the undecomposed solid. 

Early in the ablation process, we will have m, = m, = 
m = 0, until the front surface of the solid reaches its «e- 
composition temperature. Until this occurs, we will h:ve 
only one solid phase, and the heat transfer situation will be 
that of a simple semi-infinite body subjected to a sloy ly 
changing heat flux. 

For the boundary condition at the outside surface (i.c., 
Eq. [37]), we need to use relations like the ones in Equaticns 
{15-19 and 23]. In our present notation, these are 


qi(si, t) = (CHopeUe — mernr)(he — hs) — 
(0, t)]4 — AH, [11] 


oT" 1) 1... 
= = t 12 
Mor = Mm, + om 
mi = pivi(t) [14] 
¢ is the fraction of the layer material vaporized (1 — ¢ 


being eroded or otherwise lost without vaporizing). 

The transformation of Equation [41] to the doubly primed 
frame of reference (that moves with the front face of the 
layer) is accomplished by means of equations like [13 and 
14]. The result is 


oT” oT” 
[45] 
qi(si, = Bes [46] 
0) = 0 


—Ki(OT" 1/02" = qs(8s, t) + [48] 


The considerations discussed in the previous section apply 
to the solid back phase in our composite model, except for 
the boundary condition at the outside surface which is now 
determined by Equation [39] with 


qs(8s, = —K, [OT [49] 


For both phases, then, the boundary condition at the interface 
is determined by an interaction with the phase adjacent to it. 


Steady-State Solution 


In the steady state, 07'”,/Ot is zero in Equation [45] for 
the charred layer. The steady-state solution T’1«(a.) can 
then be obtained by two integrations and evaluation of the 
constants of integration from the specified boundary con- 
ditions (10). 

The first integration gives us (after combining constants) 


oT”: 
50 
|e = Bex (5 ) [50] 
bu = Ki/ po + mixCi) [51] 


which can be used in Equation [48] with the steady-state 
value of qs(ss,t) given by Equation [28]. Converted to our 
present notation, the latter is a constant given by 


oT”, pt 
ki + E —#) (v. — | + 


rr 


- 
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where (0), the ate back surface tempera- 
ture of the layer. This yields the value of the overall con- 
stant B 


B= (Comex Tix + mx SH.) exp (Le [53] 


A second integration yields the final result : : = 
(2 = A + B' exp l 51x) [54] 
where 
» 
B’ = [55] 


MN gx 


or 


AH, + SH. + C+ | 
Ds Ps 


The constant A can be evaluated at 2”, = 0 where T1*(0) = 
T,.. is the steady-state surface temperature of the layer 
A = Tux — B' [56] 
Putting x”, = Ls in Equation [54] and combining with Equa- 
tion [56], we get the relation which determines the steady- 
state wall temperature T.« in terms of the decomposition 
temperature 7)*, the ablation rates and the other model 
parameters 


Equation [57] can be solved for L* to give the result 


-1 
1 = +e] 
\ pr pi 


Tox) ps) (p./ pr 1) + Ci) 
(C.T + AH.) 


+1¢ [58] 

f 
after substitutions for B’ and B from Equations [53 and 55]. 
It is evident from this result that the dependence of the 
steady-state thickness of the charred layer on external 
conditions is expressed by the fact that it is simply inversely 
proportional to the ablation rate m.« since all other qu antities 
in Equation [58] are parameters of the model. i 

Equation [54] can therefore also be written as 


T 

[59] 
In the steady state, the thickness of the layer Lx is a con- 

stant, and we also have the equality Ss . 
(60] 


Jr, with the help of Equation [4], we can convert this to 


PsMix = PIMse = — pr)! [61] 
We also have the relation 
= Mox + Mix [62] 


The steady-state temperature distribution in the solid is 
the same as that obtained in the section “Quasi-Steady-State 
Relations” (3). In our composite model notation, this is 


T’ (2's) = Tox exp [—2' [63] 
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For this composite model, the net aerodynamic heat input 


1s 


q(net) = q" (body) + + m.AH. + (miCi + mC) X 
[7"(0) — T’(L)] + m.C.T’.(0) + mertr(he — hs) [65] 


where q¢"(body) is the rate at which heat is stored in the 
body in establishing the temperature distribution. The 
approach to the steady state can be represented by the 
condition 

lim q” (body) 


The steady-state heat of ablation is therefore 


(™*) AH, + AH. + [Te Ta) + OTe + (he — hs) [66] 


Ms* 


(he — [67] 

3* 
From the practical standpoint, it is to be hoped that 
Equations [58 and 67] can be used in calculating engineering 
estimates of important ablating properties of plastic-contain- 
ing heat shield materials. Data on the various parameters 
in these equations are presently not available. This points 
up the need for experimental studies to provide such data 
and to create a better understanding of how they are related 
to other material properties so that approximate values can 
be estimated without recourse to elaborate measurements. 
The present formulation develops the basic concepts, but is 
not necessarily complete in all details. It should serve as 
the basis for refinements of further research. 


Induction Time for the Composite Model 


The induction period, defined previously for a system 
vaporizing without layer formation, is also applicable to 
more complex models, such as the one under discussion. We 
will now deduce the corrections to Equation [34] that are 
needed for the present model of ablation. 

For this purpose, we will start with the heat conduction 
equation in a coordinate system fixed with respect to the 
ablating surface 


(ply) = + + pet) | [68 


This is a generalized equation which reduces to Equation 
[45] for the charred layer and Equation [21] for the original 
solid. In this equation, the enthalpy of the gas p,h, is 
neglected in comparison to the enthalpy p/w» of the solid 
material of the body, and v, (the velocity of the wall) is 
neglected in comparison to v, (the velocity of the gas). Con- 
servation of mass provides the relation 


Poo = (ps — [69] 


where v, is the velocity of propagation of the charred layer 
in the laboratory coordinate system. In our present model 
of the ablation process, the density of the solid body is a 
step function, equal to p, in the undecomposed material and 
to p: in the charred layer. 

It is also of interest to obtain the equation for the total 
heat H stored in the body, above that originally present. 
This quantity is related to the specific enthalpy of the body 
by the definition 


so that, by integrating Equation [68] over x, one obtains 
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The boundary condition at the ablating surface is (neglect- 
ing radiation losses) 


, oT 


(poo) no(he — hs) [72] 


= aerodynamic heat flux to sur- 
face in the absence of ablation 

amount that this heat flux is 
reduced by vaporization and 
transpiration of the solid ma- 
terial in the charred layer 

amount heat flux is reduced by 
transpiration of gaseous de- 
composition products formed 
during charring 


pwilhe + nh. h,)] 


— hs) 


It is assumed that all the ablated material is vaporized and 
that the radiation and hot wall effects are unimportant. 
Combining Equations [69, 71 and 72] gives 


vi{ pilhe = + how] pels} + 
(ps — pr)vs[no(he — hs) + how] + dH/dt [73] 


The “intercept” induction time as defined earlier and in 
Fig. 3, is 


where the integration begins at the time the heat flux q is 
turned on. The steady-state ablation velocity vx can be 
obtained from Equation [73], since in the steady-state 
dH/dt = O and v, = v; = vj. Substitution of the value of 
vx so obtained into Equation [74] gives 


| 
T= q (q vif pilh. + nlhe — hs) + how] — 
Palio + — pr) [no(he — hs) + Row]})dt [75] 


After a substitution for q from Equation [73], this reduces to 


Ifw we neglect the time dependence of the surface temperature, 
Equation [76] reduces to 


=- — pi) — hs) + hoy] Le + — [77] 


Here Lx = 


in the steady state, and Hx is the heat stored in the body in 
the steady state. Hs may be evaluated from the steady- 
state heat conduction equation (i.e., Eq. [68] with the left- 
hand side set equal to 0) 
d 
dx 


i *(v, — v,)dt is the thickness of the charred layer 


+ + (0. — ps) have =0 [78] 
or 


aT. 
+ + (ps — po)hovie] = constant [79] 


The constant may be evaluated from conditions at z = © 


and is equal to p,hyvi«; hence we have 
(poli Dix = —K, (dT x, dx) = (ps po) [SO] 


Integrating over x and using Equation [70], we get 


Vix J To 0 

=, 


Ez — + — h,) ac | [82] 


so that the induction time is 


1 T we 


Qvix To 


1 
+ (ps — pr) X 


Physically, Equation [82] shows that the induction time + 
is the sum of two terms; the first of these is simply the induc- 
tion time for a materia! which does not char; the second is 
the additional time required to supply the heat carried cif 
by the gases evolved from the charred layer. The factor 
(ps — pr)Lx in this term is the mass of the evolved gases 
per unit area; the factor [n,(h. — hs) + (1/L«) f , “(hy 

h,) dx] is the average enthalpy absorbed by these gases por 
unit mass. 

It is of interest to note that we have made no assumptions 
about the way in which the body approaches the steady 
state, and that in fact the induction time (Eq. [82]) depends 
only on the final state of the system and not at all on the way 
in which this final state is attained. This is a consequence 
of the definition of 7 which depends only on the asymptotic 
steady-state line S»(t) in Fig. 3 and not on how the asymptote 
is approached. Physically, this result is due to the fact that 
we have assumed that the amount of heat carried away 
by each bit of ablating material is independent of surrounding 
conditions, so that the total heat required to establish any 
state of the system does not depend on the path by which 
that state is reached. If radiative or hot wall effects were 
considered, or if the fraction of ablated material which 
vaporized varied with time, this would no longer be true, 
and it would not be possible to obtain as simple a formula 
as Equation [82] for the induction time. 
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Nomenclature 


steady state when used as a subscript 


Ci, = Stanton number 

C1; = heat capacity of the layer material mn g: ASeS 
in it) 

Cog = constant-pressure heat capacity of the gaseous de- 
composition products 

C, = heat capacity of the original solid 

AH. = heat of decomposition of the original solid to form 
gas and charred layer 

AH, = heat of vaporization (or decomposition) of the 
charred layer per unit mass 

€ = total emissivity of the solid surface ; 

€1 = total emissivity of the surface of the layer 

hy = specific enthalpy of the body ; 

he = stagnation specific enthalpy of gas at the edge of 
the boundary layer 

ho = specific enthalpy of the solid under initial conditions 

hs = specific enthalpy of gas at surface of body 

ky = thermal diffusivity of the layer with m, = 0 

ks = thermal diffusivity of the solid 

Ki = thermal conductivity of the layer with m, = 0 

Ks = thermal conductivity of the solid 
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‘ 
where 
i | 
@ 1 
0 q pi) (Us v1) [no(he h,)] + 
dH 
h, +—pdt [76] 
dt 
2 
| 


of 


L = thickness of charred layer ees 
Mg = mass rate of gas formation 

mi = mass rate of ablation atfrontoflayer 
Ms = mass rate of decomposition of the original solid 

Mor = Mg + om, 

ng = transpiration coefficient for gaseous decomposition 


products alone 
nT = transpiration coefficient for mixture of all gases 
produced (i.e., from decomposition of solid and 


vaporization of layer) 


* = steady-state heat of ablation 
q(z, t) = heat flux in stationary frame of reference at station 
x and time ¢ 
gs(8s, b) = heat flux into undecomposed solid in the stationary 
frame of reference 
q's(“'s, ) = heat flux in primed frame of reference at station 


x’; in the solid at time ¢ 
q’"1, ) = heat flux in doubly primed frame of reference at 
station x”; in the layer and at time ¢ 


7 


face of the original solid) 

T’;(x's, t) = temperature in primed frame of reference at sta- 
tion x’; in the solid at time ¢ 

T",,X", = temperature-position coordinate axes moving with 
the surface of the layer 

T” \(x",, t) = temperature in doubly primed frame of reference 


at station x”, in layer at time ¢ 


Ue = velocity of gas at edge of boundary layer 

Ug = velocity of the ablating gas 

vy = velocity of front surface of layer (velocity of the 

wall) 

Vs = velocity of back surface of layer (front surface of 
original solid) 

x1 = coordinate along surface normal in layer in sta- 
tionary frame of reference 

ent = coordinate along surface normal in 7”;, X’; frame 
of reference 


a = coordinate along surface normal in 7"), X”; frame 


of reference 


Po = density of the body (equal to ps in original solid 
| and p; in charred layer) 
pl = bulk density of charred layer sn, 
ps = density of original solid References 
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of Heat Transfer ‘and 
Friction Coefficients for Helium 
Flowing in a Tube at Surface 
sabe — up to 5900 R 


HE IMPORTANCE of nuclear reactors as a power 

source for aircraft and space vehicles has stimulated 
interest in convective heat transfer from high temperature 
surfaces. The resulting large ratio of fuel element tempera- 
ture to working fluid temperature means a large variation 
in the properties of the fluid, which influence the heat transfer 
characteristics. Some work has been done with wall tem- 
peratures up to 3460 R using air (1). Heat transfer co- 
efficients for helium flowing through a carbon tube with a 
maximum inside surface temperature of 5040 R and a corre- 
sponding gas temperature of 4640 R are presented in (2). 

In order to extend the range of surface temperature re- 
ported in (1) and the ratio of surface to fluid bulk tempera- 
ture of (2), the present investigation was undertaken. In this 
investigation the wall temperatures were increased to the 
limit of the tungsten and molybdenum test sections, and the 
ratio of surface to fluid bulk temperature was kept as high as 
possible. 


Experimental Apparatus 


Arrangement 


A schematic diagram of the arrangement of the test section 
and equipment used in the investigation is shown in Fig. 1 
The mixing tanks and test section were housed in a steel 
containment tank, which was evacuated to about 50 u of 
mercury. This was done to minimize oxidation of the re- 
fractory metals of the test sections. Pressures lower than 50 
x of mercury were possible but would have resulted in evapora- 
tion of the hot metals. 
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of average heat and friction coeffic local heat coefficienis 
were made with helium flowing through electrically heated smooth tubes with length-diametcr 
ratios of 60 and 92 for the following range of conditions: 
4533 R, Reynolds number from 3230 to 60,000, heat flux up to 583,200 Btu per hr per ft? of heat 
transfer area, and exit Mach number up to 1.0. The results indicate that, in the turbulent ranze 
of Reynolds number, good correlation of the local heat transfer coefficients is obtained when the 
physical properties and density of helium are evaluated at the surface temperature. 
heat transfer coefficients are best correlated on the basis that the coefficient varies with [I 
(L/D)~:"] and that the physical properties and density are evaluated at the surface temperature. 
The average friction coefficients for the tests with no heat addition are in complete agreement with 
the Karman-Nikuradse line. The average friction coefficients for heat a addition are in poor agre e- 


ment with the accepted line. 


Average surface temperature from 1457 ‘o 


The average 


Test Secti ; las 
Test Sections 

Since the object of this taswtigtion was to obtain hext 
transfer data at as high a surface temperature as possible, 
this experiment embraced a number of materials problems. 
There are only five known metals with melting points above 
5000 R: Molybdenum, osmium, tantalum, rhenium and 
tungsten, in order of increasing melting points. All five 
materials are subject to rapid oxidation when heated in air. 
For reasons of availability and high melting point, molyb- 
denum and tungsten were chosen as test section materials. 
Since molybdenum tubes were available commercially, they 
were used for the first high temperature tests. Several tubes 
were made by disintegrating holes through tungsten rods, 
which were then ground to a uniform wall thickness. Nickel 
entrance and molybdenum exit flanges were used to connect 
the test sections to the mixing tanks. The molybdenum test 
section was equipped with a bellmouth entrance. A right- 
angle edge approach was used for the tungsten test section 
entrance. In addition to the refractory metal test sections, a 
well-instrumented Inconel test section was used to make 
several runs to check the high temperature heat transfer 
results. 

The dimensions of the various tubes used as test sections 
in this investigation are shown in Table 1. 


| 


Table 1 Dimensions of tubes used as test sections in 
this investigation 


Heat transfer 


Tube ID, OD, length 
material in. in. L, in. L/D 
molybdenum 0.191 0.250 11.5 60 | 
| tungsten 0.125 0.250 11.5 92 | 
Inconel 0.124 0.188 11.25 91 H 
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Fig. 1 Schematic diagram showing arrangement of test 
apparatus 


J 
Instrumentation 


‘he outside wall temperatures of the molybdenum and 
tungsten test sections were measured with platinum-platinum- 
13-per cent-rhodium thermocouples spotwelded along the 
length of the test section and located as shown in Fig. 2. For 
test section temperatures beyond the range of the thermo- 
couples, a small target optical pyrometer was used. The use 
of the optical pyrometer is discussed in the Appendix. Some 
experimentation with tungsten-molybdenum thermocouples 
was conducted, but the extreme brittleness of the tungsten 
wire after welding made the use of this thermocouple im- 
practical. The Inconel test section wall temperatures were 
measured with Chromel-Alumel thermocouples. 

The gas temperature was measured before entering and 
after leaving the test section with platinum-platinum- 
rhodium thermocouples located downstream of the baffles 
in the two mixing tanks. 

The radiation shield and outside surface of the Transite 
insulating cylinder were instrumented with platinum- 
platinum-rhodium thermocouples. Molybdenum wires were 
spotwelded along the length of the test section to measure 
voltage drop as a function of distance from the entrance. 
Static pressure taps were located in the entrance mixing tank 
and in the entrance and exit flanges of the molybdenum test 
section. Voltage and static pressure tap locations are shown 
in Fig. 2. The flanges of the tungsten test section were not 
instrumented with pressure taps because of the difficulty in 
getting small holes through the molybdenum flange and the 
walls of the tungsten tube. 

The molybdenum test section was used in obtaining both 
friction factor and local and average heat transfer data. 
Since the tungsten test section was not instrumented with 
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‘Fig. 2. Typical test section. Thermocouple, optical pyrometer, 


__- voltage tap and pressure tap locations shown are for molybdenum 
test section. (All dimensions in inches) 
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Fig. 3 Variation of thermal conductivity and absolute viscosity 


of helium with temperature 
Method of er 


Helium Properties 


The variable physical properties of helium used in calculat- 
ing the Nusselt, Reynolds and Prandtl numbers are shown in 
Fig. 3 as a function of temperature. In the absence of 
experimental data, the curves in Fig. 3 were calculated by 
the methods shown in (3, 4 and 5). The value of specific 
heat cp is 1.24 Btu/(lb)(R) and is constant; the ratio of 
specific heats y is 1.667, and the gas constant RB is 386 
ft-lb/(Ib)(R). 


voltage or pressure taps, it was used only in obtaining average a _ 
heat transfer data. The range of conditions for which data - 
were obtained is summarized in Table 2.  & ‘The average friction coefficient f was calculated by the 
Table 2 Summary of range of conditions for which data were obtained ; 
Bulk Reynolds Average surface 
Test section number temperature Exit Mach flux, 
material Re» number -Btu/(hr) (ft?) 
t molybdenum 3750-60 , 000 530 0.09-0.97 adiabatic 
molybdenum 3230-25, 230 1457-2334 0.21-0.95 73, 100-300 , 000 
tungsten 5375-11 ,000 1905-4533 0.72-1.0 169 , 600-383 , 500 
Inconel 15,340-35, 100 1606-2057 eer 272 ,400—583 , 200 
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method used in ( 1) wherein 


4(L/D)(pavV?/2g) 2(L/D)G? 


where the density pa» was evaluated at the average static 
pressure and temperature of the gas as follows 


1 (pi t+ 

The average friction coefficient was also calculated with the 
density evaluated at the film temperature, as shown in (1). 


Heat Transfer Coefficient 


The average heat transfer coefficient ha, was computed 
from the experimental data by the relation 


Q wep (Ts — 


~ §7.—T,) S(T. — T) [3] 


hav 


where the temperature of the gas 7, was taken as the arith- 
metric mean of the total temperatures at the entrance 7; 
and the exit 7. of the test section. The average surface 
temperature 7, was taken as an integrated average of the 
local outside wall temperature minus the temperature drop 
through the wall. The temperature drop through the wall 
was calculated by the following equation, which is derived 
in (6) 


where 


= average outside wall temperature a 
thermal conductivity of the tube material 
inside and outside radii of the test section, respec- 
tively 
In this investigation the temperature drop calculated by 
Equation [4] was very small compared with the difference 
between the average inside wall temperature and the helium 
bulk temperature. 


«>, 


Because the test section temperatures were very high, i 
was expected that the heat loss to the surroundings would be 
large. For most runs, less than 20 per cent of the heat gencr- 
ated was lost, according to calculations by the following 
equation 


where 


Q: = heat loss to the surroundings 
Q. = electrical heat input 
Q = heat transferred to the gas 


Most of the heat loss could be accounted for. 

It was possible to calculate local heat transfer coefficie::ts 
by evaluating the various heat losses and then making a he:t 
balance for each increment. The heat loss at the ends of t':e 
test section and the radial heat loss for each increment w: re 
calculated. These losses accounted for more than 80 per cent 
of the total heat loss calculated by Equation [5]. Each 
local heat loss and the end losses were increased by the ratio 
of total heat loss to the sum of local and end heat loss:s. 
The rates of incremental electrical heat generation and heat 
conduction into and away from each increment were calc ii- 
lated and used along with the heat losses to set up a heat 
balance for each increment. From this heat balance the 
rate of heat transfer to the gas could be calculated for each 
increment. Starting at the entrance of the test section, the 
bulk temperature of the gas leaving the first increment could 
be calculated. This was repeated for each succeeding incre- 


‘ment. The local bulk temperature and local surface tempera- 


ture were used to calculate a local heat transfer coefficient. 
The local and average heat transfer coefficients were used to 
caleulate Nusselt numbers. 


res 
Results and Discussion | 
Axial Wall Temperature Distribution 


Three representative axial outside wall temperature dis- 
tributions for the tungsten test section showing thermocouple 
and optical pyrometer temperature measurements appear in 
Fig. 4. The experimental data as indicated are summarized 
in Table 3. (Continued on page 889) 


Q-/S, Q/S, w, 
Run Btu/(hr) (ft?) Btu/(hr) (ft?) lb/hr 


2374 
2976 
4305 
5349 
3202 
5674 


~ 


93, 200 
132,200 
187,000 
266 , 000 
191,000 
362,500 


73, 100 3. 
109, 200 
173, 500 14. 
208 , 000 19 
159, 200 7 
300,000 5.9 


to orien 


A 


227 ,000 
461, 500 
753 , 000 
1,145,000 


169, 600 
307 ,000 
383, 500 
362, 500 


324,000 
320, 500 
520,000 
522,000 
600 ,000 
599 , 000 
652,000 


275,400 
272,400 
457, 200 
486 , 800 
551, 700 
516 , 900 
583 , 200 


Table 3 Experimental results 


Pi, 
lb/ft? abs 
Molybdenum test section 


Tungsten test section 


Inconel test section 


P2, 
lb/ft? abs 


1040 
1155 
1464 
1732 


| 
| 
R R R R 
ag 099 580 1514 1047 2259 
ace pele 2 056 565 1130 848 1738 
3 2063 560 1021 791 1497 
ae bate 4 4 2088 560 978 769 1457 
TA. 5 2093 572 1383 978 2257 
6 2198 560 1193 871 2334 
615 1465 1905 
if 7.2 620 1690 2499 j 
7.9 650 2278 3418 
1 4.6 745 —-2720 4533 
a 11 530 1230 880 1606 
12 8.5 535 1310 923 1688 
ne 13 16.8 530 1192 861 1652 
ro 14 18.4 525 1165 845 1584 
15 17.3 ve 525 —«-1300 912 —-1869 
16 15.3 530 1350 940 1921 
17 15.7 526 1430 978 2057 
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